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Abstract 



We study frustrated, two-dimensional, quantum antiferromagnets in the 
vicinity of a quantum transition from a non-collinear, magnetically-ordered 
ground state to a quantum disordered phase. The general scaling properties 
of this transition are described. A detailed study of a particular field-theoretic 
model of the transition, with bosonic spin-1/2 spinon fields, is presented. Ex- 
plicit universal scaling forms for a variety of observables are obtained and the 
results are compared with numerical data on the spin-1/2 triangular antifer- 
romagnet. Universal properties of an alternative field-theory, with confined 
spinous, are also briefly noted. 
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I. INTRODUCTION 



Three kinds of frustrated 2d 
mentally and theoretically. First 
VCh, VBr2, CqEu, NaTi02 etc [|] 
1973 when Anderson and Fazekas I 



There has been a remarkable recent revival of interest in the low-energy properties 
of two-dimensional {2d) frustrated quantum antiferromagnets. In part, this interest was 
triggered by the discovery of strong magnetic fluctuations in the high-T^ superconductors; 
however, frustrated magnetic systems are interesting in their own right, in the light of nu- 
merous theoretical predictions on the nature of disordered ground states in quantum spin 
systems 

systems have been studied intensively, both experi- 
are antiferromagnets on a triangular lattice such as 
Theoretical studies of such antiferromagnets go back to 
J] first suggested that ioi S = 1/2, quantum fluctuations 
may be strong enough to destroy the classical 120° ordering of Heisenberg spins. Though 
most of the subsequent numerical and analytical studies do indicate [[7| the presence of long- 
range order at zero temperature (T), these studies also show that quantum fluctuations 
are quite strong. 

A second frustrated system is the antiferromagnet on a kagome lattice. It is believed 
to describe the second layer of ^He on graphite |T0| and SrCr^^xGai^j^xOig and related 
rT|. The effects of quantum fluctuations in kagome antiferromagnets are far 

and numerical studies of S* = 1/2 systems support a 
quantum-disordered ground state at T = |T3|,3 . Besides, large S kagome antiferromagnets 
display the Villain order from disorder phenomenon : in the semiclassical approximation, 
they possess a strong 'accidental' degeneracy which is lifted only by the zero-point motion 
of quantum spins [p!2|,p!5| . Tunneling between a sequence of nearly degenerate ground states 



compounds 
stronger than in triangular ones [112 



(which differ in energy only due to quantum fluctuations), may also contribute substantially 
to the reduction in the strength of the large 5* long-range order 
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Finally, there are also studies of antiferromagnets on the square lattice which are frus- 
trated by adding second and third neighbor couplings PJT7[] . These systems show interesting 
phases with incommensurate, planar, spiral correlations. 

A key feature of the local spin correlations in the systems above, which will be crucial in 
our analysis, is that they are non-collinear. Unlike the unfrustrated square lattice, the spins 
are not locally either parallel or anti-parallel to one another. The analysis in this paper will 
mostly assume that the spins are coplanar although this second restriction is mostly in the 
interests of simplicity. 

So far we have discussed the situation at T = 0. Experiments, however, are performed at 
flnite T when thermal fluctuations are also present. The effects of thermal fluctuations for 
2d Heisenberg systems are well known [|T^ - they destroy long-range magnetic correlations 
at arbitrary small T. Suppose, flrst, that the ground state is nearly perfectly ordered. It is 
clear, then, that at small T, thermal fluctuations will be signiflcantly more important than 
quantum fluctuations, and the low-T behavior will be predominantly classical - the primary 
effect of quantum fluctuations will be a renormalization of the couplings at T = 0. This 
is the low-T "renormalized-classical" 
later observed 



regime which was studied in detail in Ref p5|, and 



^ in a number of experiments on undoped square-lattice antiferromagnets 
at sufficiently low T. Consider, next, the physics when the system is quantum-disordered at 
T = 0. Then all thermally induced fluctuations are suppressed by a (presumed) spin- gap at 
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low enough temperatures and and the low-energy dynamics is purely quantum mechanical - 
this is the "quantum-disordered" regime. 

However, there is a third, intriguing possibility which arises when the ground state of 
the system is not too far from a T = 0, second-order quantum transition between the 
magnetically-ordered and quantum disordered states. Then it is easily possible to find the 
so-called "quantum-critical" regime where classical and thermal fluctuations are equally im- 
portant. This is a high temperature regime with respect to any energy-scale which measures 
the deviation of the ground state of the antiferromagnet from the quantum transition point; 
on the magnetically-ordered side a convenient choice for this energy-scale is a spin stiffness, 
Ps- However it is also a low temperature regime with respect to a microscopic, short- distance 
energy scale like a nearest-neighbor exchange constant, J. If the couplings are precisely crit- 
ical, then the quantum-critical region stretches down to lowest T - this is unlikely to be 
realized in antiferromagnets without fine-tuning of an external parameter e.g. pressure or 
doping. However, even if the system is not precisely at the critical point, but T is larger 
than Ps on the ordered side, or a corresponding energy scale A on the disordered side, we 
still observe essentially quantum-critical behavior because at such T we effectively probe 
the system at scales where it does not know on which side of the transition it will end up 
in its ground state. However, if the long-range order at T = is very well established (or, 
if on the quantum-disordered side, A is very large) the condition ksT > ps {ksT > A) for 
quantum-criticality may conflict or interfere with ksT < J and the the quantum-critical 
behavior can be overshadowed by nonuniversal short-range fluctuations. Thus we require 
that Ps (A) be reasonably small, and then then we may expect to observe quantum-critical 
behavior at T smaller than J. 

In a recent publication with J. Ye [^, two of us considered whether a quantum-critical 
region exists in the square-lattice S = 1/2 antiferromagnet. We computed various exper- 
imentally measurable quantities such as the uniform susceptibility, the correlation length, 
the dynamic structure factor, and the spin-lattice relaxation rate for antiferromagnets with 
coUinear spin correlations. We found reasonable agreement between the quantum-critical 
results and the experimental data [20| on La2-xSrxCu04^ and with numerical results on 
S = 1/2 antiferromagnets [^]. We argued, therefore, that this system is quantum-critical 
at intermediate temperatures. In frustrated 2d systems, quantum fluctuations are likely to 
be far stronger. It is therefore reasonable to expect that quantum-critical behavior may be 
observed in frustrated systems as well. In the present paper, we will present detailed predic- 
tions about the quantum-critical properties of frustrated antiferromagnets with non-collinear 
correlations, to help elucidate this possibility. 

The study of quantum-critical behavior is not the only purpose of our analysis. We 
will also consider the behavior of various observables in the renormalized-classical region. 



Previous studies in this this region were performed by Azaria et. al. [|18|, who focused 
on a renormalization group analysis for the correlation length. Below, we present, for 
the first time, expressions for the uniform susceptibility and dynamic structure factor of 
renormalized-classical, non-collinear antiferromagnets. 

An important issue, which makes a study of non-collinear antiferromagnets considerably 
more difficult than collinear ones, is that the nature of the quantum-disordered phase and the 
universality of the transition are not well established. The large- Sp{N) theories have 
argued that the quantum-disordered phase of non-collinear antiferromagnets has deconfined, 
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spin-1/2, bosonic spinons. In this paper, we derive a macroscopic field-theoretical model 
which has the same behavior, and study the universality class of the transition between 
a quantum-disordered phase with deconfined spinons and the magnetically ordered state. 
We find, quite generally, that such a transition is in the universality class of the 0(4)- 
vector model in spacetime dimension D = ?>. This result agrees with the semiclassical 
renormalization group analysis of the magnetically ordered side m D = 2 + e dimensions of 
Azaria et. al |T3]. We will then go on to determine numerous universal, finite temperature 



properties of such antiferromagnets. These properties have many striking differences from 
those of the coUinear antiferromagnets [^] which possessed confined spinons. Note however, 
there are other treatments of the transition in non-collinear antiferromagnets which do 
not have 0(4) exponents at D = 3. We will review these in Appendix A and show that they 
in fact have confined spinons. The universal magnetic properties of these approaches differ 
only in a minor way from those of Ref and will therefore not be discussed in any detail. 

We will begin in Section by defining carefully, and with considerable generality, the 
order parameter of coplanar antiferromagnets p3|j26| , p!8[| . We will also express the staggered 
dynamic susceptibility in terms of correlations of the order parameter. We will continue our 
general discussion in Section [1B| where we will present universal scaling forms for nearly- 
critical coplanar antiferromagnets. These scaling forms follow from not much more than 
the presence of hyperscaling and a dynamic critical exponent z = 1. On the magnetically- 
ordered side, the entire dynamic staggered and uniform susceptibilities will be argued to be 
fully universal functions of five parameters characterizing the ground state: Nq, the order- 
parameter condensate, the two stiffnesses py, and the two susceptibilities xy, X± (defined 
more precisely below). Similar results will hold also on the quantum disordered side. We 
emphasize that none of the results of these two sections make any specific assumptions on 
the universality class of the transition. 

In Sections ^^^^ present explicit computations of the universal scaling functions 

using a particular (we think likely) field-theoretic model of the transition. This approach has 
deconfined spin-1/2 spinon excitations in the quantum disordered phase, which lead to many 
interesting observable consequences. Section VII will compare some of the above results 
with available numerical results for the S* = 1/2 triangular Heisenberg antiferromagnet; this 
comparison will use some new results on the 1/5* expansion of this model which are obtained 
in Appendix B. 

Our main conclusions will be reiterated in Section [VIII] . The contents of Appendix A 
were noted above, and some technical details will be presented in Appendix C. 



A. Order parameter and other observables 

For simplicity, we will restrict our discussion to antiferromagnets with Hamiltonians of 
the following form: 

H = Y.J,,Si-S, (1.1) 

where the Sj are spin 5* operators on the sites i,j of a regular two-dimensional lattice, 
and the Jij are the exchange integrals. The Jij respect the symmetries of the lattice, and 
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are short-ranged, although not necessarily nearest neighbor. The strength of the quantum 
fluctuations will depend on the value of S and on the ratios of the Jij, and will determine 
whether the ground state is magnetically ordered or quantum disordered. 

In the following, it will be convenient to think of the Sj not as quantum operators, but 
as spacetime-dependent fields in a path integral over imaginary time r. We will restrict our 
analysis to antiferromagnets in which the strongest fluctuations are well described by the 
following hydrodynamic parametrization 

Sj(r) = ni(xj, r) cos(2Q ■ xj + n2(xi, r) sin(2Q ■ xj (1.2) 

where rii, 112 vary slowly on the scale of a lattice spacing, but are always orthogonal: rii ■ n2 = 
for all Xj, T. The ordering wavevector 2Q may be commensurate or incommensurate with 
reciprocal lattice vectors, but must not be such that ( |1.2D makes all the Sj collinear with 
each other. Thus the square lattice with 2Q = (vr, 7r)/a is excluded (a is the nearest neighbor 
spacing), as is any ferromagnetic state (Q = 0). The triangular lattice and certain kagome 
lattice antiferromagnets with 2Q = (Svr/S, 87r/v^)/a, or square lattice antiferromagnets with 
incommensurate Q are however included. Kagome antiferromagnets with more complicated 
local correlations, which are nevertheless coplanar, will also be described by our universal 
results, but are not considered explicitly for simplicity. The parametrization (|1.2| ) also 
implies that the spin orientations are always locally coplanar. In fact, even antiferromagnets 
with non-coplanar correlations can be analyzed by a straightforward extension (not described 
here) of our results. The key restriction is that the correlations are non-coUinear: we however 
assume coplanarity for simplicity. 

As is well-known p3| , p^ , |26[] , we can identify the pair of vectors rii, 112 as the order 
parameter of the antiferromagnet; below we will discuss an equivalent complex matrix order 
parameter, Qa,i3, which is computationally somewhat more convenient. On the magnetically 
ordered there will be a spin-condensate which, we assume, satisfies 

N'o = (ni)|=o = (n2)ko (1-3) 

Our analysis will rely heavily on a spinor parametrization of the vectors rii, 112. This is 
most directly introduced by the Schwinger boson representation of the spin operators 

Sa = \hial,hp (1.4) 

where a = x,y,z, a, P = 1, 2, and the a" are the Pauli matrices; site and time dependence 
of the fields is implicit. It turns out that the hydrodynamic form (|1.2|) is related to the 
following parametrization of the b: 



baiir) = (2,(x„r)e^Q-» +2£,^;2;(x„r)e-*Q-''') (1.5) 

where e is the antisymmetric tensor and the Za are slowly varying complex fields satisfying 
the following normalization at some scale A: 

N 

Y^\za\' = N (1.6) 

a=l 
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with N = 2 (we have introduced the variable N in anticipation of the generahzation below to 
arbitrary A^). The renormalization factor Zs accounts for the fluctuations at scales shorter 
than A. Inserting ( p. .51 ) in ( |1.4| ) and comparing with ( |1.2| ) we obtain 

1 • a n ^7^ 

n2a + mia = -^Sa^Z^Cr^pZp (1.7) 



It is easy to check that this satisfies ni ■ n2 = and (|1.3| ). Notice that order parameters 



fields are quadratic in z, this is consistent with the identification of the z quanta as 5 = 1/2 
bosonic spinous. The composite character of the order parameter was also noticed (for 
= 2) in Ref. pf. 



Some key properties of the above parametrization deserve notice. First, the question 



of gauge invariance. As is well-known, the Schwinger boson decomposition (|1.4| ) demands 
that the physics be invariant under the U{1) gauge transformation h — > 6**^6. However the 
continuum parametrization (|1 . 5| ) 'breaks' this gauge symmetry [Q. Alternatively stated, if 
the z fields are slowly varying in one particular choice of gauge for the 6, they will have 
forbidden rapid variations for most other gauges. Thus, simply by focusing on a long- 
wavelength theory of the z, we have 'broken' the gauge symmetry. There is however, a 
remnant Z2 gauge symmetry that must be kept track of: notice that the transformation 

2;(x, r) ^ ?7(x, r)2(x, r) : 77 = ±1 (1.8) 

leaves all the spin operators invariant. All observables must be invariant under this Z2 
gauge transformation. All of these features are consistent with earlier large A^ theories of 
frustrated antiferromagnets which found breaking of U{1) gauge invariance down to Z2 
in all non-coUinear antiferromagnets. 

Consider, next, the symmetries any effective action for the z must satisfy. It must clearly 
be invariant under any global SU{2) spin rotation z Uz, where U is an SU{2) matrix. 
More interesting, however, is the behavior under lattice translations x — > x + y. The 



spin-rotation invariance of Ti, and the parametrization ( |1.5| ) are consistent with this only if 
the action is invariant under the global transformation 

z e-'^-^z (1.9) 

where y is any near-neighbor vector. For the triangular lattice, this demands that the action 



be invariant under the Z^ symmetry ||2^ z exp(±z27r/3)z, while for incommensurate spiral 
states it is effectively equivalent to a global U{1) symmetry. In practice we will find that the 
consequences of the Z^ symmetry are essentially identical to the larger U{1) symmetry, and 
we will therefore simply refer to this lattice symmetry as a U{1) symmetry. It is important, 
however, not to confuse this global, lattice, U{1) symmetry, with the U{1) gauge symmetry 
discussed above. Thus the effective action for the z field should possess a global SU (2) x U{1) 



symmetry [29|; for general A^ this will be a SU{N) x U{1) symmetry. 

An important observable which will characterize the antiferromagnet, is the staggered 
dynamic susceptibility Xs defined by 



Xsik, iujn)5ab = Ylj ^'^ {Sia{r)Sjb{0)) exp [-i ((k + 2Q) ■ (x^ - xj) - u;„r)] 

(1.10) 
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at the small momentum k away from 2Q and Matsubara frequency ujn- The sums over i,j 
extend over all the Ns sites of the system, and Vg is the volume per spin (e.g., Vs = a/3/2 
for triangular antiferromagnet). 

The physically measurable retarded staggered susceptibility can of course be obtained 
by the usual analytic continuation to real frequencies. Inserting (|1.4| ) and (|1.5|) in ( |1.10D , we 
find 

1 ^ r rh/ksT 

Xsik^iu^n) = ^^^^^^^ Y.J d'x rfr(g„^(x,r)g;^(0,0))e-('^-™) (1.11) 

where N = 2, and the symmetric order-parameter Qap = is given by 

Qaf3 = ^J^ZaZp. (1.12) 

Note that it has A^(A^ + l)/2 different complex components, and is invariant under the Z2 
gauge transformation (|1.8|) . It transforms under SU (N) x f/(l) 2 Young tableau under 

SU{N) and as charge 2 under U{1). Again we have introduced an A^-dependent notation 
to facilitate the generalization to arbitrary N. The equation ( |1.3| ) for the magnitude of the 
order parameter can also be expressed in the general form 

^0 = E \{Qo.p)t=o\ (1-13) 

a(3=l 

We also quote for reference the relationship, special to iV = 2, between the tensor order 
parameter Qap and the vectors ni, n2, which can be deduced from (|1.7|) and ( |1.12|) : 

Q ^ }^ f -n2x + in2y - irii^ - riiy + iui^ \ (I U) 

2\ n2z + iriiz n2x + in2y + iui^ - ^ly / 

We will find it convenient to express many of our results in terms of the dynamic, stag- 
gered, structure factor which is the Fourier transform of the spin-spin correlation function 
in real time t 

/roo 
(fx / dt(5z(x,t)5^(0,0))exp-i(kx-cut) (1.15) 

This is of course related to the staggered susceptibility defined above by 

2h 

In addition to the order parameter, the uniform magnetization density, M(x), is an 
important hydrodynamic variable. Its fluctuations decay slowly due to the conservation law 
for the total magnetization. It is defined by 

M(x„t) = ^S,(t), (1.17) 

where gfis/^ is the hydromagnetic ratio. Its diffusion is measured by the uniform spin 
susceptibility defined by 

^ r roc 

Xu{k,uj) 6ab = -j-J dt{[Ma{^,t),Mb{0,0)])exp-i{k^-ujt) (1.18) 
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B. Scaling forms 



We will now consider the properties of the non-collinear antiferromagnets in the vicinity 
of a second-order quantum phase transition from a magnetically ordered to a quantum disor- 
dered ground state. We will try to keep the discussion in this section as general as possible, 
independent of any specific field theory for the transition. The results of this subsection 
will follow from some fairly general scaling assumptions, rather similar to those applied to 
collinear antiferromagnets in Ref . A primary assumption will be that the quantum tran- 
sition has dynamic critical exponent z = 1. Explicit computations of the scaling functions 
and exponents will be performed in the subsequent sections using a particular deconfined- 
spinon field-theory of the transition. A confined-spinon field-theory will be briefly considered 
in Appendix A; its properties are also consistent with the scaling ideas of this section. 

Let us assume that the T = transition occurs as some coupling constant g is varied 
through a critical value g = gc, and the magnetically ordered state occurs for g < g^. 

We present first the scaling properties for g < g^. We expect that the condensate Nq will 
vanish as 

No-^igc-gf (1.19) 

where /9 is a universal critical exponent. A second characterization of the ordered ground 
state is provided by the spin stiffnesses p\\ and p±: these measure the energy cost of twists in 
the plane and perpendicular to the plane of the order parameter, respectively. In the presence 
of hyperscaling (which we assume), we expect that both these stiffnesses will vanish as 

P±,P\\--igc-gr (1-20) 

where the u is the usual correlation length exponent (this formula is special to two dimen- 
sions). Further, the ratio of these two stiffnesses will obey 

lim = (1.21) 

gygc p± 

where Tp is a universal number. In a similar manner we can consider the two uniform mag- 
netic susceptibilities X± defining the response of the antiferromagnet with infinitesimal 
anisotropy to uniform magnetic fields perpendicular and parallel to the plane of the order 
parameter, respectively (note the inversion in the order of 'parallel' and 'perpendicular' !). 
In a 2; = 1 theory their scaling properties are identical to those of the spin stiffnesses, and 
possess an associated universal ratio T^. The subsequent sections of this paper consider a 
field theory in which Tp = T,^ = 1 exactly; in Appendix A we briefly consider a model with 
different universal ratios. In all cases it is useful to consider the dimensionless numbers j/p, 

Vx 

^^^Pl^r^ . = (1.22) 

P± x± 
which measure the deviation of the stiffnesses and susceptibilities from the universal ratio 



at the critical point; clearly yp,yx — as 5^ — > (^c- Finally, as in Ref |25|, we also need 
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the following dimensionless ratios which measure the wavevector, frequency, and stiffness in 
units of the absolute temperature 



J- hc±k _ hu NksT 
ksT ksT 4:7ip± 

The numerical factor of An is for future notational convenience, and the spin-wave velocities 
c±, C|| are of course given by c]_ = p±/x± and c| = P\\/x\\- The factor in xi has been 
inserted because oc in the large limit, and so ensures that xi remains of order unity 
in this limit. 

Now, following arguments closely related to those in Ref . , we may conclude that the 



response functions of nearly-critical antiferromagnets obey the following universal scaling 
forms 

X.(*.^) = ^ ^) [^) (1.24) 

Xu{k, u) = ^ ^bT $!„ (k, uJ, xi, Vp, (1.25) 

Here and Si are completely universal functions of their dimensionless arguments 

and there are no non-universal scale factors anywhere. The exponent f] is related to the 
order parameter exponent f3 by the hyperscaling relation 

2^={1 + T])u. (1.27) 

From the above scaling relation and (|1.19|) and (|1.2CI|) we see that the prefactors of all 
the scaling functions remain finite all the way up-to g = Qc, or xi = oo. Further, all 
scaling functions are defined such that they remain finite as xi ^ oo when we will also 
find i/p^y^ —>■ 0. The universal functions $1^ and Si are related by the fluctuation-dissipation 
theorem Si = Im$is. As in the argument xi determines whether the system is better 
described at large scales by a quantum-critical {xi ^ 1) or a renormalized-classical {xi -C 1) 
theory. 

Strictly speaking, the leading scaling properties of the observables are obtained at 
yp — Vx ~ ^1 because these ratios are associated with irrelevant operators. However many 
long-distance properties are sensitive to the precise values of the spin-stiffnesses and sus- 
ceptibilities. Thus these operators are actually dangerously irrelevant, and it necessary to 
consider many observables as full functions of yp and y^ 

Parallel arguments can be applied to the quantum disordered state with g > g^. We 
assume that this state has low-lying quasiparticle excitations with non-zero spin, character- 
ized by an energy scale A, which propagate with a velocity c. In the model considered in 
the subsequent sections we will have spin-1/2, bosonic quasiparticles above a gap A; there 
are however other possibilities, one of which is discussed in Appendix A. We expect that A 
will obey 
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A - QcT 



(1.28) 



near the critical point. We also introduce the dimensionless ratio 

-2 = ^ (1.29) 

which is the analog of the xi on the ordered side. There is now no need to consider the analogs 
of the Up, as these will be truly irrelevant (as opposed to dangerously irrelevant) on the 
disordered side. The observables of the nearly-critical, quantum-disordered antiferromagnet 
obey 

XsiK 1^^^ (^^^ $2. (f, cU, X2) (1.30) 

Xu{k,Uj) = (yt) {k,L0,X2j (1.31) 



S{k^.) = nA(^] (^) (1.32) 



he ( knT^ ^ 



Again $2s, ^2u and S2 are completely universal functions. The pref actor A is related to 
quasiparticle amplitude(s) and vanishes as 

Ar^ig-g^r (1-33) 

The precise definition of A requires a normalization condition on $2s which will be discussed 
later. 

Before closing this section, we briefly introduce the scaling functions of some other im- 
portant observables which can be deduced from the ones above. We restrict ourselves to the 
ordered side; the extension to the disordered side is straightforward. The scaling function 
for the spin correlation length is 

C' = ^^Xix,,yp,y^) (1.34) 

The static uniform spin susceptibility at g < Qc behaves as 

Xu{T) = {^y^BT 0(xi,y„yj (1.35) 

The local structure factor Sl{u) is given by Sl{u!) = J cPk S{k,uj)/A7r'^. The contribution of 
Xu to Sl{uj) is subdominant and Sl is given simply by a momentum integral of the staggered 
susceptibility. This integral is always ultraviolet convergent (because the intermediate states 
in S{k,u!) are all on-shell) and is dominated by k less than about 1; we have therefore 

where Ki = J (fk Ei/Ati^. The small frequency limit of S^iuj) is directly related to the 
spin-lattice relaxation rate 1/Ti oc Sl{P — > 0). We will also discuss static structure factor 
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S{k) = J duj S{k,Lj) /27r. The frequency integral is divergent at the upper cutoff if fj > 1, 
whence S{k) is non-universal - this will be the case in our model. 

In the subsequent sections we will obtain explicit expressions for the scaling functions 
introduced above in the renormalized-classical and quantum-critical regions. We will use 
a new deconfined spinon field theory which will be introduced in Section along with a 
expansion which will facilitate our computations. Properties of the quantum-disordered 
phase will also be discussed. 



II. EFFECTIVE FIELD THEORY: DECONFINED SPINONS 

The main aim of the remainder of this paper is to illustrate the general scaling ideas 
discussed above in the framework of a specific field theoretical model of the quantum tran- 
sition. An important property of field-theory we use is that results of the spacetime 
D = A — e, and D = 2 + e, expansions on it are all consistent with each other; we will 
consider only the 1/N expansion here. 

A significant reason behind the choice of our particular model is that it possesses decon- 
fined spin- 1/2 spinon excitations in the quantum disordered state. This is then consistent 
with the S'p(A^)-large N prediction of Ref. on non-collinear antiferromagnets. Further, 
our approach will allow us to explore some of the observable consequences of these novel 
excitations. 

We begin with some discussion on the role of the Z2 gauge symmetry of (|1.8| ). The 
crucial role of this gauge symmetry was noted in was emphasized to us at an early stage 
by N. Read and was also noted in Ref [0]. Our main assumption will be that the Z2 
gauge symmetry can be entirely neglected in the continuum field theory. In other words, 
configurations with a non-zero local Z2 flux remain gapful across the transition. The Z2 
gauge fluxes are in fact present in the cores of vortex lines (in spacetime) associated with 
homotopy group tti{S0{3)) = Z2 of the true 5*0 (3) order parameter [^. We assume that 



these vortices remain confined across the transition and that the Z2 gauge charge of the 
z field is globally defined (the 2;-field configuration around a vortex is double- valued) . 
Under these circumstances we may simply write down a continuum Landau- Ginzburg field 
theory for the 2;-field. Implicitly, this procedure implies that we are not distinguishing 
between SU{2) and S0{?>) symmetries. 

We will now write down the most general action consistent with the SU{N) x f/(l) 
symmetry discussed before. Rather than using a soft-spin Landau- Ginzburg approach, we 
find it more convenient to use hard spins satisfying ( p..6| ); this modification is however not 
crucial and completely equivalent results can be obtained by the former approach. To 
second-order in spatial gradients this yields the following effective action 



S = I cPxdr — 

fJ.=X,T 



df^Z^d^Za {z^d^Za d^^Z^ Za) 



(2.1) 



where a = 1 . . . N, and gx, Qt, 1x, It are coupling constants. Any of these couplings can be 
varied to tune through the quantum transition - we will use 

g = 9x- (2.2) 
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Simple considerations presented in Section |T| below show that these coupling constants are 
given by 

-JL -JL _PlzA _4zA fo^^ 

and and x° are the bare values of two spin stiffnesses and spin susceptibilities. For 
simplicity, throughout the paper we define transverse and longitudinal susceptibility without 
a factor g^s/fi'- A more detailed consideration of the values of p, x is presented in the next 
section. 

The effective action S can also be explicitly derived from microscopic considerations. 
Using the continuum parametrization in (|1.5|) it is not difficult to show that the long- 



distance limit of the Sp{N) theories of Refs and |jT2[ is described precisely by S. The 
same parametrization can also be used on the semiclassical approach of to obtain S. 
Finally, we also explicitly derived the effective action for N = 2 from the general 



macroscopic approach of Ref [23 



Some critical properties of S can be immediately deduced. By a simple power-counting 
argument in D = 4 — e dimensions it can shown that the 7^ couplings are irrelevant at the 
critical point. An identical result can also be obtained hj D = 2 + e analysis parallel to that 



of Ref |]T8[. We will also explicitly show the irrelevancy of the 7^ in the 1/N expansion below. 
(None of these arguments of course exclude the possibility that a large bare value of 7^ may 
have more fundamental effects. In fact, at 7^ = —N, S actually becomes the U{1) gauge 
invariant CP^~^ model, which is then a model for quantum phase transitions in collinear 
antiferromagnets. We will not consider the possibility of these large 7^ complications in this 
paper.) 

It is therefore useful to begin the analysis by considering 5 at 7^ = 0. It is easy to verify 
that now S has its internal symmetry enlarged from SU{N) x U{1) to 0{2N). Further, the 
spacetime theory is Lorentz invariant. Finally, this theory has py = p± (and similarly for x) 
and so we have Tp = = 1 exactly 

The exponents appearing in the scaling functions are now all properties of the well-known 
0{2N) fixed point, and we quote for reference to order (see also Appendix C) 

Note that the exponents f] and (3 are associated with the composite field Qaf3 and thus differ 
from the usual rj, (3 for vector fields. Thus fj is quite close to unity at large A^, while the 
corresponding r] which appears in collinear antiferromagnets is almost zero. 

Let us now consider how the 7^ variables break the Lorentz and 0{2N) symmetry. Te- 
dious but straightforward computations show that the terms proportional to the 7^ transform 
under a single irreducible representation of 0{2N) - the one labeled by a Young tableau of 
2 rows and 2 columns. It is therefore not necessary to decompose the 0{2N) structure of 
the operator. However the 7^ terms are irreducible under the Lorentz group - there are the 
spin-0 and spin-2 pieces 

71 = (27, + 7,)/3 

72 = (7x-7r)/3. (2.5) 
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The terms associated with the 71 and 72 are now completely irreducible under 0{2N) and 
Lorentz group, and will therefore have their independent crossover exponents, 0i and 02 
respectively, measuring their irrelevancy. In other words, near the quantum fixed point, the 
fully renormalized spin stiffnesses and susceptibilities should obey 



Pll - P± 

P± 
X\\ - X± 

x± 



7i(6)-'^^+72(e 



J) 



7i(0)"^^- 272(0)-'^^ (2.6) 



where O is the Josephson correlation length measured in lattice units. To leading order in 
7 we also have from ( p.6|) for the spin-wave velocity difference. 



c± 2 



72(0)"^^ (2.7) 



As g,, approaches g"^, O behaves as O ~ (l-fl'^/S'^i) " (clearly, g^/ g^ = gr/gr)- In section|Tl 
we will find the following expansion result for the renormalization-group eigenvalues 
attracting the 71^2 to the fixed point 

32 112 



III. CONSERVED CHARGES AND CURRENTS 

This section will present the computation of the spin stiffnesses and uniform spin suscep- 
tibilities both at T = and in the quantum critical region of the deconfined spinon action 
S. The calculation will be carried out to order We will show how one can obtain 

renormalized stiffnesses in the ground state by doing calculations in the symmetric phase 
at T ^ 0. A computation of stiffnesses directly in the ordered phase is performed in the 
Appendix C. 

The stiffnesses and uniform susceptibilities are all response functions associated with the 
conserved charges and currents of S. We will therefore begin by studying the SU{N) x 
U{1) symmetry of S. The conserved charges and currents can be determined by the usual 
procedure of evaluating the change in the action under an infinitesimally small symmetry 
transformation with a spacetime dependent angle. The results are conveniently expressed 
in terms of the A^^ — 1 traceless Hermitian SU (N) generators which we choose to satisfy 

Trace (T'^T'') = ^5'^'' (3.1) 

Then the currents associated with the SU (N) symmetry can be written as 

K; = -- {z^T-d.z - d.z'^T-z) - ^ [z^d^z - d,zh) [z^T^z] (3.2) 

The index /i extends over x, r and is, strictly speaking, a conserved charge density - in 
this section we will use the term 'current' to generically refer to both charges and currents. 
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We will also not explicitly display the spacetime-dependence of the fields. The current 
associated with the U{1) symmetry is 



9n 



(3.3) 



Our intention is to express the fully renormalized stiffnesses in terms of SU {N) and U (1) 
current-current correlators, and so we need the appropriate Kubo formula. To derive this 
formula it is convenient to introduce vector potentials which linearly couple to the conserved 
currents above, and examine the response of the system to these vector potentials. Let us 
consider first a SU (N) vector potential A^^. This modifies the action to 



S' ^ — I d^xdr \ f . 
91^-' L ^ 



'd^ + tAlT-]z 



- (df,z^ - iAlz^T"") z) ' 



It is then not difficult to obtain the response of the free energy F = — log / Vze ^' 
external vector potential. Doing the algebra we find 



- + — {z^ (T'^T^ + r''T«) z)^ + ^ Iz^T'zz^T^ 



(3.4) 



to the 



(3.5) 



Again space-time dependences have been suppressed, and the two fields inside the correlator 
are at different spacetime points. A very similar analysis can be carried out for a [/(I) vector 
potential a^, and we find 



5^F 
Sal 



2A^(1 + 7m) 
9n 



(3.6) 



Now we change tracks and evaluate the response of the system to these vector potentials 
in an entirely different way. Let us assume that we are on the ordered side {g < gc), and are 
able to integrate out all the fluctuations, including the amplitude fluctuation modes in the 
direction of the condensate. We then obtain a fully renormalized action for the spin-wave 
fluctuations. Let this effective action have the following form 



(3.7) 



Here Z is a A'^-component complex vector of unit length which yields the local orientation 
of the condensate. The factor of 2 in F is introduced for further convenience. Let the 
condensate point in some fixed direction Zq = (1, 0, 0, 0, . . .). We now look at small variations 
about this direction as in 



Z ^ Zo + {ia, TTi 17^2, TTs + iT^i, ■ ■ 



(3.8) 



the factor of 1/2 is present because Z is a spinor and rotates by only half the angle of the 
observable order parameter. The effective action for this variation is 
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dPxdi 



2N-2 



1=1 



(3.9) 



By the definition of the stiffnesses we identify pi^, P2x as the two spin stiffnesses of the spin 
wave modes: 

Plx = P± , P2x = P\\ (3.10) 

Also the stiffness to twists in the time direction gives us the uniform spin susceptibility: 

piT = X± , P2T = X\\ (3.11) 

Now let us look at the response of F to the presence of an external SU{N) vector potential, 
while the condensate is non-zero. Doing the same analysis as before we obtain 



6F = 2 (fxdr 



p,,AlAlzlT'^T'Z, + (p2M - Pi,) A^A"^ [ZlT'^Z,] [ZlT'Z, 



7]rpb I 



(3.12) 



For a fixed condensate Zq this result will depend upon the orientation of the SU{N) rotation 
A^. However if we place the system in a box of large, but finite length in the p direction 
the response of F is clearly proportional to 5"* because no symmetry can be broken (for the 
case p = T this equivalent to having a small finite temperature T oc L~^). Thus we should 
replace each T°-T^ factor by its average over all the generators of SU{N) - it is crucial that 
we average over all the generators, and not over different orientations of the condensate. We 
will then need the identities 



Ar2 



- — y zIt^T^Zo = — 

-l4" ° 2N 



2N{N+1] 



(3.13) 



which are actually true for any complex unit vector Zq. These identities can be easily 
established by considering explicit forms for the T'*. So finally, combining ( |3.12|) and ( p. 131) , 
we can determine the response of F to the SU{N) vector potentials at an infinitesimal 
temperature: 



S'F ^ , 2_ Nx± + XII . S'F 

N {N + l) ' 



,ab 2 Np^ + p\\ 



5"' — 



N {N + l) 



(3.14) 



We will evaluate the left-hand side using (|3.5|) and thence obtain an expression for the above 
linear combination of the stiffnesses. 

We still need a second linear combination - this is of course provided by the U{1) currents. 
An exactly parallel computation can be done for the response to the U{1) vector potential 
a-r - in this case we find 



6al 



4X11 



5^F 
Sal 



(3.15) 
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Combined with ( |3.6| ), (|3.5| ) and ( p.6| ) we now have reduced determination of the spin stiff- 
nesses and susceptibihties to evaluation of the correlators in ( p.5|) , (|3^) at T = 0. This 
calculation will be carried out in Section [ill A| to order 1/N . 



These methods can also be used to obtain the temperature dependence of the uniform 
spin susceptibility Xu{T). By analysis similar to that in Ref ||25| it is not difficult to show 
that 



where there is no summation over a. This computation will be considered in Section |TTB. 

We conclude with a note on the nature of the 1/A^ expansion of S. We found that a 
properly renormalized theory for the scaling functions can only be defined if we restrict with 
the leading terms in an expansion in powers of 7^: we shall therefore do a double expansion 
in powers of 1/A^ and 7^. This expansion is most easily done be treating the effects of 7^ 
perturbatively - i. e. without introducing a Hubbard - Stratonovich decoupling of the quartic 
term. 



A. Spin-stifFnesses and susceptibilities at T = 



Below we will need the form for the vertex function associated with the anisotropic term 
in the action. In the momentum space we have 



^ a,l3{kl,^^ ^3,/i; 



AN 



(3.17) 



where a and /3 number the components of the 2;-field. The diagrammatic representation for 
the current-current correlation functions is shown in Fig.|l|. At = 00, one can neglect 
self-energy and vertex correction within a bubble; however the renormalization due to F 
generally cannot be neglected because the summation over the components of the z-field in 
the extra bubble associated with F yields a factor of A^ which cancels out the 1/A^ factor in 
( p.lTp . However, a simple inspection of the diagrams shows that the effects of F are relevant 
at A^ = CX3 only for the U{1) correlator, while for the SU{N) currents, the side vertices in 
the bubble contain sign-oscillating T matrices, and the summation over the components of 
z— field gives only a factor 0{1). As a result, we find at A^ — *■ 00 and in the limit T — 



2N 



6^F 



1 







(3.18) 



and 



6^F 



2iV l-T 



1+7,. 



(3.19) 



where g"^ = Qc and = c\gc- Clearly from ( p.l4| , |3.15|) , the r.h.s. in ( |3.18|) and (|3.19|) 
are also the values of 4p_|_^ and ipy respectively. Note that, as one might expect, only py 
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acquires a correction due to 7, while p±_ remains the same as in the isotropic case. Also note 
that ( p. 19 ) is indeed consistent with ( p.6|) and establishes that 0i = 02 = 1 at = 00. 

We now describe the 1/A^ corrections. Obviously, we have to consider the self-energy and 
vertex corrections within a bubble in Fig.|I|, and the renormalization of the vertex function 
r itself. The latter however is again relevant only for U{1) response, while for SU{N) 
response, the leading effect of 7^ is itself of order 'Jfi/N and there is no need to consider 
the renormalization of F to order 1/A^. The computation of the SU{N) response therefore 
requires less efforts, and evaluating the diagrams in Fig.|l| with F given by ( |3.17|) , we obtain 



2N 



6^F 



2N6' 



ab 



l+7M/(2iV) 



9n 




1 - 



5V 
9u. 



2u 



(3.20) 



where (yfj^ 



1 + 7^/2 A^) and u = 1 — 16/37r^A^ is the critical exponent for the correlation 

length. 

Our next step will be to calculate, with logarithmic accuracy, the renormalized value of 
F as T — >■ 0. We will then use the result to compute the U{1) response to order 1/A^. The 
diagrams which contribute to the vertex renormalization to order 1/A^ are shown in Fig. |^. 
The internal part of each diagram contains two Green functions and the polarization operator 
- this combination produces logarithms after integration over intermediate momentum and 



frequency in 2+1 dimensions |2^. The evaluation of diagrams is tedious but straightforward, 
and after doing the algebra we obtained that the momentum dependence of the vertex 
remains the same as in ( p.l7|) but 7^ changes to 7^'^-^ where 



128 

15^ 
48 



log(l -^a;M) +7 



32 



16 



log(i - 9x1 gc 



48 , ,^ , A 16 , , 



(3.21) 



Substituting the renormalized vertex into the bubble diagram for t/(l) response and per- 
forming also self-energy and vertex renormalizations within each bubble in the way described 



in we obtain to order 1/A^ 



6al 




+ 7: 




(3.22) 



Eqns ( ^.20 ) and ( |3.22| ) can now be combined to obtain transverse stiffness to order 1/A^ 



4pj 



2N 
9tM 



1 + 



2nJ 



1 - 



5V 



-ye// 

2N 



1 - 



5Vi 



(3.23) 



Finally, using (|3.22| ) and ( |3.23| ), we obtain the result for (py^^ — p±,^)/p±,^ to order 1/A^. 



Reexpressing 7^ in terms of correct renormalization group invariants 71 and 72 of (|2.5|) , and 
exponentiating logarithmic terms, we obtain the 1/A^ results for the crossover exponents (pi 
and 02 that were given in ( p.8|) . Our value of 02 coincides with the result by Lang and 
Ruhl who computed anomalous dimensions of tensor fields of arbitrary rank for critical 
0{2N) sigma models. On the other hand, there do not seem to be any other computations 

of 01. 
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B. Uniform susceptibility 



The calculation of the uniform susceptibility at small but finite T and arbitrary p± is 
essentially the same as that of SU{N) response at T — >■ 0; only the summation over frequency 
should not be substituted by the integration. Doing the same calculations as have let us to 
( p.20| ) but at finite T, we obtain to first order in 

Xu{T) = (^)' m (l + ^ m) (3.24) 

Here 



x{T) = x{T = 0) + 6x{T) (3.25) 

where 



We expect that higher-order corrections to (|3.24| ) will only change 7t- to ^^J^ . Notice that 
at T — i> 0, we recover a result consistent with ( p. 151) and ( |3.16|) : 

The temperature dependent piece SxiT) in ( |3.24[ ) is precisely 1/2 of that in the isotropic 
0{2N) sigma model with iV— dependent spin- wave velocity c*. At = cxd we have from ||2^ 



5x(T) = {ksT /2ticy) /(xi), where numerically f{xi) is close to 1 for all kBT/pj_. We will 
describe the structure of 1/iV corrections to 6x and the value of c* in the following sections: 
the 1/A^ results are of a rather different physical form depending upon whether ksT ^ p± 
or ksT <^ pj_. We will therefore consider the expressions for Xu{T) and other observables 
separately in the renormalized-classical and quantum-critical regions. 
We now begin our discussion of various low-T regions. 



IV. RENORMALIZED-CLASSICAL REGION 



This section will present expressions for different scaling functions in the renormalized- 
classical region, fc^T ^ p±. Under this condition, the low-temperature behavior is related 
to the low-energy fluctuations of the macroscopic order parameter of the ground state and 
is therefore essentially classical. Indeed, this is true only for fluctuations at sufficiently large 
scales when typical energies huj ~ hc±k <^ ksT, and one need consider only the oj^ = 
term in the summation over Matsubara frequencies. At larger k, quantum fiuctuations are 
important, and at k > C,J^, antiferro magnet possesses D = 2 + 1 critical spin fiuctuations 
We will consider this critical behavior in the next section, and here focus on the situation 
at small hc±k < ksT. As in unfrustrated antiferromagnets, there are two different low-T 
regimes already in the classical region, because the actual (thermal) correlation length ^ is 
exponentially large when ksT <^ p^, and one can have either k^ <^ 1 oi k^ ^ 1 P5[ . 
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Physically, the crossover at fc^ ~ 1 is between the regime where the ordering is destroyed 
by classical fluctuations and the dynamics is purely relaxational {kC, < 1), and the regime 
where classical fluctuations are weakly damped propagating gapless spin-waves {kC, > 1). 
Below we will see how the spin structure factor changes in passing from one regime to the 
other. But flrst we consider the behavior of the correlation length. 



A. Correlation length 

As in the collinear case, we deflne the correlation length from the equal-time, long- 
distance, exp(— r/^) decay of the spin-spin correlation function. From our previous dis- 
cussion, especially from and (|1.12|) , it is clear that the Fourier transform of the 



spin correlator is related to the polarization operator rather than to the Green function 
of the z~ fleld. At = oo, spinous behave as free particles, and their propagator is 
Go{k,iLj) = l/i^k"^ + uj'^ + ml), where mo is the mass of the 2;— field, which at = 00 
coincides with the inverse correlation length of the 0{2N) model. We then obtain 

g(r)oc / , ,,,,, ;"5f,f' , oce--- (4.1) 



[(q + k/2)2 + m2] [((q-k/2)2 + m, 



We see that in this limit, the actual correlation length, ^, is precisely l/2mo. We now proceed 
to finite A^. To first order in we have to consider self-energy and vertex corrections 

within a polarization bubble. A simple inspection of terms shows that while self- 
energy corrections renormalize the spinor Green function, and hence ^, vertex corrections 
only modify the overall factor in the correlation function and do not affect the exponent in 
the decay rate. In other words, to first order in 1/N, the actual correlation length is still 
precisely a half of that for the fields, and we therefore only have to locate the pole in 
the zero-frequency part of the field propagator. For the isotropic case (7^ = 0), such 
calculations have already been performed in . Here we have to consider also the effect of 



the 7^ anisotropy. It is not difficult to check that the anisotropic term contributes to the self- 
energy to first order in and therefore affects at this order the constraint equation which 
in essence is the equation for ^. The 7-dependent self-energy piece can easily be calculated 
because the only nonvanishing contribution to order comes from the diagram in Fig.|I|. 
We obtain 

J:,{k,^u;) = ^u^ + ^clk' (4.2) 



where Cq = c° = \J p\/x\. Let us first keep only anisotropic self-energy term. Substitut- 
ing ( [4.2|) into the constraint equation ( |1.6|) and performing the momentum and frequency 
summation, we obtain 



where m is the full mass for the z— field, c is the linear combination of the two spin-wave 
velocities which we will compute below, and is the same as in (|3.2CI| ). We now observe that 
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the r.h.s. of ( [4.3|) can in fact be reexpressed in terms of the fully renormalized transverse 
and longitudinal spin- stiffness. Using ( p. 20 ) and ( 3.23|) , we find 



log 



2n 



hem 



2p± 
N 



1 



1 



(pii - p±y 



2{N + 1) 



(4.4) 



Our next step is to determine how ([4.4]) is modified by other corrections. We first 
consider the change in the r.h.s. of ( [4.4|) as T ^ 0. At 7^ = 0, earlier calculations to order 
p5[ have shown that the only temperature- independent modification of the constraint 

(N 



equation is the renormalization of the coupling constant gx in (|4.3| ) to gx (N — l)/N. This 
renormalization can effectively be regarded as the wavefunction renormalization of the z- 
field, such that each z— field propagator acquires a factor Z = {N — l)/N. Physically, 
this renormalization is related to the fact that the solution of the constraint equation at 
arbitrary small T and finite p exists only for > 1, while for = 1 (i.e., for the XY 
case), a single gapless spin- wave mode has no partner to interact with. Consider now the 
7-dependent piece in the r.h.s. of ( ^^) . Clearly, it should also acquire an extra factor 
similar to the renormalization og gx- It is difficult however to determine explicitly the 1/A^ 
renormalization of the anisotropic term because the anisotropic vertex itself has a factor of 
1/A^. On the other hand, the form of the wavefunction renormalization seems quite plausible 
from a physical perspective, and we assume below, without proof, that it remains the same in 
the anisotropic case as well. Simple considerations then show that 'jx should be substituted 
by ^x/Z = 'yxN/{N — 1). We then obtain, keeping only temperature- independent corrections 
in the r.h.s. of 



log 



2ps 



27r 



hem N — 1 



(4.5) 



where = Ps/x-, ^"^^ Ps and x are given by 




N {pw-PlY 



2{N^ - 1) ; 

N {X\\-X±y 
2(Ar2 - 1) XL , 



(4.6) 



Finally, we collect all temperature-dependent 1/A^ corrections to ( [4.4|) using the same proce- 
dure as for the 0{2N) model p5[. These corrections include double logarithms in the form 
ksT log log fcsT/m, and regular O^ksT) terms. Double-logarithms eventually give rise to 
the temperature-dependent prefactor in ^. Assembling all contributions, we finally obtain 
for the actual correlation length in frustrated antiferromagnet 



1 ^ he 

2 ^k^ 



\N-i)kBr 

47rp, 



l/2(Af-l) 



exp 



4:7iPs 



{N-l)kBT 



(4.7) 



where 36 



l/2(Af-l) 



X r(l + l/2(A^-l)) 



(4.8) 
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We see therefore that, to first order in 7, the expression for the correlation length is the same, 
up to a factor of 1/2, as in the 0{2N) isotropic sigma-model with effective spin-stiffness Aps 
and spin-wave velocity c. The factor of 4 in ps merely reflects the difference between the 
definitions of the coupling constant g in (|2.1| ) and in the 0{2N) a— model. At the same 
time, the overall factor of 1/2 is a signature of deconfined spinous. For the physical case of 
= 2 we have 



Ps 



3PX + -Pii 



(4.9) 



To first order in 7 we also have c = 2c_l/3 + C||/3. The T dependence in ( [4.7|) then agrees with 
the two-loop renormalization group calculation of ^ performed by Azaria et al [13]. They 



also obtained the two-loop expression for the correlation length in a classical model, valid 
at arbitrary ratio of the two bare stiffnesses, and argued that the result for the quantum 
case at arbitrary 7^ will be the same if expressed in terms of the fully renormalized p± and 
P||. Our analysis shows that this universal behavior of the correlation length certainly exists 
to first order in 7^, but we have no proof that the universality persists at arbitrary 7^. In 
any event, the analysis presented here is valid close to the critical point when higher order 
corrections due to anisotropy are small. 



B. Uniform susceptibility 

The result for Xu{T) valid at arbitrary ratio of T/p± is given by ( p.24 ). We now use the 



results of Ref 25 and obtain 



N-1 keT 
N 2^ 



(4.10) 



where c is given by ([4.6|) . We expect that this result will hold at arbitrary A^. For the 
physical case of A^ = 2 we then obtain using ( [4. 101) and (|3.24| - |3.26|) 




x± + 



Attc^ 



(4.11) 



In Sec |VI1] we will apply our result for Xu{T) to the S* = 1/2 Heisenberg antiferromagnet on 
a triangular lattice. 



C. Staggered susceptibility and structure factor 



The key input for this subsection is our observation, in Eqns. ( p..7|) , ( p..l2| ), that the 
hydrodynamic order-parameter variable for frustrated antiferromagnets is a composite op- 
erator made of two z— fields, and spin-spin correlation function is related to the polarization 
operator of spinous. At A^ = 00, we use (|1.11|) and the results of Appendix C, and express 
Xs{k,iuj) as 



Xs{k,iu) 



4pi 



n(fc, iu) 



(4.12) 
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Here A'^o and p\_ are the fully renormalized values of the on-site magnetization and spin- 
stiffness at T = 0, and 11 is the polarization operator which for hck, hu < ksT is given 
by W 



n(/c, iuo) 



k^T log [{k^ + 0^ + ^(P + tu2)2 + 4Pm2) /2kmo 



TT 



(4.13) 



As before, mo is the mass of the 2— field at A^ = 00 which to this order is also the inverse 
correlation length for the z— field, and oj = oj/cq. We see that at small k < mo, liik^ioj) ~ 
ksT / Aurn^ and hence Xs{k-, ioj) oc T^^, where ^ is the actual correlation length. At the same 
time, at fc^ ^ 1, the logarithm in the numerator of ( [4.13| ) cancels the overall factor of T, 
and we obtain Xs{k,iuj) oc 1/fc^ as it should be in the Goldstone regime. 

We now consider how this simple behavior is modified by 1/A^ corrections. A simple 
inspection of the 1/A^ terms shows that the diagrams which contribute to the renormalization 
of Xs are the same as in Sec iV A - they include isotropic self-energy and vertex corrections 
within a polarization bubble, and also corrections due to the 7^. Let us first consider the 
isotropic case. The self-energy corrections to the z— field at 7^ = were studied in p5[ . 
They give rise to a renormalization of the mass and of the bare stiffness, and also yield an 
overall thermal renormalization factor Afe for each Green function. For /c ~ cj/c^ 3> m, this 
renormalization factor is 



A. 



A^- 1 
A^ 



1/2 



loglksT / {hc±m) 
log[V^^ + fn'^/ IT"' 



-1/2{N-1) 



(4.14) 



At k 



~ m. 



the logarithm in the denominator is a number of the order of one, and we have 



Xk = [kBTiN-l)/A7rp^^ 



l/2(Af-l) 



[1 + 0{1/N)). Further, it is not difficult to check that the 



vertex renormalization within a bubble also gives rise to logarithmic terms. Evaluating the 
corresponding diagram in Fig. ^, to accuracy 0{1) and exponentiating the result, we obtain 
another renormalization factor (k, which at A; ^ m and to order 1/A^ is simply = Af. 
Collecting both contributions, we then obtain 



Xs{k,iuj) 



Nl 
4pi 



At U{k,tu) 



(4.15) 



Finally, we consider the effect of the anisotropic term to first order in 1/A^. Clearly, there 
are self-energy corrections to the z— field propagators which eventually change p± to ps and 
c± to c. Besides, the anisotropic terms give rise to vertex corrections. We didn't perform 
actual calculations of the vertex corrections, but on general grounds it is likely that, to 
order 0{1/N), they can be absorbed into the renormalization of A"o. We then assemble 
(PD, (|1|), ( CT) and (|415| ) and obtain 



Xs{k, 



lUO] 



PsiN-1) 



kBT{N - I) 
Aups 



-, (iV+l)/(7V-l) 



e f{ki,uji/c) 



(4.16) 



where the overall factor is chosen such that /(O, 0) = 1. It follows from ( [4.16| ) that at finite 
A^, Xs(0,0) oc T^N+i)/{N-i)_ r^YAs result is likely to be vahd at arbitrary A^. For the 
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physical case of = 2, it reduces to xTi^^ 0) °^ T^C"^ - this is substantially smaller than the 
naive mean-field result Xs(0, 0) oc T^^. 

The behavior of f{x,y) at intermediate x,y = 0(1) is rather complicated, chiefiy 
because the spin-wave velocity also acquires a substantial downturn renormalization at 

= 0(1) However at fc^ ~ cj^/c ^ 1, the velocity renormalization is irrelevant 

and using ( [4.14|) and (|4.15| ) we obtain 



^'"•'''= (fTi) ^ <5 log(.Vy=))'"«""'-" (4.17) 



Substituting this result into ( |4.16| ), and using the fact that at ^ 1, logx ~ log^, we 
obtain to first order in 1/A^ 

„(M-)^(^)fj^ (4.18) 

We now demonstrate that at any N, this expression is nothing but the rotationally-averaged 
spin-wave result for the ordered SU{N) x U{1) antiferromagnet. Indeed, using ( |1 . 1 1| ) , ( |1.12| ) 



and ( p.8|) we find that the total number of transverse spin waves in the ordered state is 
Nsw = 2A^. That ( [4.1^ ) is the averaged spin-wave result now follows from ( |1 . 1 1| ) and 



the fact that each transverse spin-wave mode at T = contributes a spin-wave factor 
NQ/2ps{kP' + uj"^ /c^) to Xs (see Appendix C). For the physical case of = 2, the averaging 
factor is Nsyj/N{N + 1) = 2/3, as it should be. 

For experimental comparisons, it is also useful to have an expression for the dynamical 
structure factor S{k,u) defined in ( [1.15| ), and static structure factor 

Sik)= I^Sik,co) (4.19) 

As before, we will be interested in the behavior of S{k, u) at scales much larger than 
the Josephson correlation length. At such k, quantum fiuctuations are irrelevant, and using 
( p..l6| ) and (|1.26| ), we can conveniently reexpress S{k,uj) as 



where Si is related in a straightforward manner to the universal function Hi introduced 
earlier in ( [L.26|) . Below, we will restrict consideration of S(/c,C(j) to the frequency range 



C(j ~ c/^, which is relevant for experimental comparisons in the renormalized classical region. 

Consider first, the asymptotic behavior of S(A;,Ci;) at large momentum k^ 3> 1. For such 
k, we found above that 1/A^ corrections to the polarization operator are not singular. For 
a qualitative analysis, we can then safely restrict ourselves to = cx), when the imaginary 
part of the polarization operator is given by 

lmU{k,u) = ^ [ [(1+n,, (5(e,, +e,, -M- (^(e,, + M) + 

{nq, - Uq^) (5(6^2 - eg, - huj) - 6{eq^ - eg, + fiuj))] (4.21) 
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Here Uq is a Bose function and = {hc)^q^ + ^ ^- At ck ^ cj, the only contribution 
to Im n(/c, a;) comes from the second piece in ( 4.21j ), which describes coUisionless Landau 



damping. Doing the integration, we obtain a simple result 

Note that as expected, 'E.i{k^uj) scales linearly with uj. 

We turn now to a discussion of smaller k. The corrections to ( 4.21|) include the terms 



similar to A^, C,k above, which grow logarithmically with decreasing k and eventually change 
the temperature dependence of Hi at k^ ~ 1. Moreover, at such momenta, the damping of 
excitations becomes comparable to the real part of the quasiparticle energy, and we cannot 
simply restrict ourselves to coUisionless Landau damping. We did not perform explicit 
calculations of Si at intermediate fc, but for an estimate we can rely on the results 
of Ref p7| , |25|] for the momentum dependence of the damping of excitations in the 0{2N) 



sigma-model. Combining these results with the expressions ( [4.16| , [4.17| ) for the real part of 
Xs, we obtain 



2/{N-l) 



2i(A:,a.)oc-p^p (4.23) 

Here •jkuj is the damping of z-field excitations given by [H'il] 7fcw oc 
hckk {ksT / p^gY log p'^/ksT, and the momentum-dependent spin-stiffness and spin-wave 
velocity are introduced as another way to account for the logarithmical terms in ( [4.17[) : 



, {N - l)knT 



Ps = log k^; {cuY cc pt (4.24) 

(We note in passing that at A;^ ^ 1, we have with the logarithmical accuracy p^ = ps,Ck = c.) 
At k^ = 0{1), we have p^ oc T, Ck oc \/T at arbitrary N, and hence, our final result 

For = 2, we have Ei{k, u) oc u T^'^. 

Finally, we consider the static structure factor, S{k). A simple inspection shows that the 
frequency integral in ( |4.19| ) has two basic contributions. One comes from large uj where the 



system is -D = 2 + 1 critical, while the second comes from huj < ksT where the fluctuations 
are essentially classical. We will see in the next section that at large uj, ^.{k^uj) behaves as 
l/tu^"*' where fj is given by ( p.4| ). We found earlier that f/ > 1 (at least, at large A^, and hence 
the frequency integral over quantum fluctuations explicitly depends on the upper cutoff in 
the theory. We will discuss nonuniversality in S{k) in more detail in the next section. In 
the renormalized-classical region however, the correlation length is exponentially large and 
we may expect that the dominant contribution to S{k), which scales as still comes from 
the frequences uj oc where fluctuations are essentially classical. For such frequences, the 
rescaling factor between Im x{k,uj) and S{k,uj) is 2/(1 — ^^/^bT-^ ^ 2kBT/hcj, and we have 
simply S{k) = kBTxs{k,0), where Xs{k,0) is given by ( |4.16| ). At = we then obtain 
8(0) oc T2^/(^-i) e^. For A^ = 2, this reduces to 8(0) oc T'^^. 
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D. Local susceptibility and spin-lattice relaxation rate 



Another experimentally measured quantity is the momentum-integrated dynamical struc- 
ture factor S{uj) = J d?k S{k,uj) / Att"^ . Unlike S{k), this observable is universal in 2D as can 
easily be seen from ( [4.22|) . It is also not difficult to show that for uj ~ g^"^, the integration 
over momentum is confined to /c ~ where we can use the estimate ( [4.20| , [4.25| ) for 
S{k,Lj). We then obtain 

Further, the — >• limits of S{k,uj) and are related to the transverse (1/Ti) 

and longitudinal (I/T2) relaxation rates for nuclear spins coupled to electronic spins in the 
antiferromagnet. We have 

1 , r d^k .2 



Ito yj-^AfS(l,.) (4.27) 
2''I"Jt^Y fll^^txlik..) (4.28) 



T2J NMR "^^0 V/IC/ J A-K^fi 

where and A^. are the hyperfine coupling constants (with the dimension of energy). They 
generally tend to some finite values as A; — > 0. The factors of 2 appear because fiuctuation 
modes near Q and —Q equally contribute to relaxation rates. The temperature dependence 
of 1/Ti then immediately follows from the result (|4.26|) for S{uj). For = 2 we obtain 



1 ^ (Ay fksTV" ^^^^^ 



Ti \h J c \ ps 

An exactly parallel computation can be done for the spin-echo decay rate I/T2, and the 
result is (for general A^) 




(4.30) 



For A^ = 2 this yields T^^ oc T^i. 



V. QUANTUM- CRITICAL REGION 

We now consider the results for the quantum-critical region where 47rps < NksT. Under 
this condition, the relevant scale for fiuctuations is given by T itself and both quantum 
and classical fiuctuations are equally important (i.e., at relevant energies, Bose functions 
are 0(1)). Our first observation in this region concerns the role of the anisotropic {i.e. 7^ 
dependent) terms in the action. The scaling hypothesis predicts that any scaling function 
near the quantum transition should depend on the dimensionless ratio ^j/L^ where is 
the Josephson correlation length, and = hc/ksT is a finite length in the imaginary 
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time direction at g = Qc- We have shown above in Sec |TT| that at T = 0, anisotropic 
corrections had a form "^fii^j / ci)~'^^''^ where both crossover exponents are clearly positive 
and even larger than 1 at finite N (see Eqn. ( p.8| )). We therefore expect that the leading 
anisotropic corrections deep in the quantum-critical region will scale as ^^{kBTa/hc)'^^'^ with 
positive 01^2, i-e. they will be subdominant at low T compared to the leading terms in the 
scaling functions. Clearly then, the quantum-critical behavior will be the same as in the 
isotropic 0{2N) sigma-model. The anisotropic term in the action will however renormalize 
the spin-wave stiffnesses and velocities in the subleading terms in the full scaling functions, 
which describe deviations from the pure critical behavior. These terms will be calculated 
in this Section only at = oo, at which order the anisotropic term in the action does not 
contribute. We will then assume, without proof, that the renormalization due to 7^ leads 
to the same effective ps and x given by ( [4.6| ) as the renormalized-classical expressions. On 
general grounds, this is likely to be the case because the corrections to the pure quantum- 
critical formulas account for the crossover to the renormalized-classical region. However, as 
we said, explicit calculation of the subleading terms at finite N has not been performed. 

We emphasize that even in the absence of the anisotropy, the scaling properties of spin 
correlators are quite different from those for unfrustrated antiferromagnets simply because 
each spin component is a bilinear product of the 2;— fields. We now consider separately the 
behavior of various observables. 



A. Correlation length 

The expression for the correlation length follows directly from the observation that the 
spin propagator is a convolution of two Green functions for z-fields. An analysis, similar to 
that for the renormalized-classical region, shows that vertex corrections in the polarization 
bubble do not effect the form of the exponential decay of correlations, and therefore the 
actual correlation length is again exactly 1/2 of that for the 0{2N) sigma-model. Specifically, 
we obtain 

^T) = I ^ X,{oc) [1 - Kx^'/-' + . . .] (5.1) 

where we defined Xi = NksT / Aixps-, and v is the exponent for the Josephson correlation 
length given by (p.4|) . The values of Xi(oo) and k were found earlier [^: Ai(oo) = 
6(1 + 0.1187/A), where 6 = 2 log (1 + v^)/2 = 0.962424, and k = 2/v^ + 0{1/N). 



B. Uniform susceptibility 



We continue with the response to the uniform magnetic field. As in the renormalized- 
classical region, we use the general result (p.24|), but now the temperature dependent piece 



in X is dominant, and to order 1/N the universal function for Xu defined in (|1.35|) is given 
by 



\/5 , v^+1 
— log 

27r ^ 2 



1 - 



0.31 

1^ 



+ axi + 



(5.2) 
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where a = 0.8 + 0{1/N). This is indeed a half of the susceptibihty for 0{2N) square- 
lattice antiferromagnet. At = 2, we obtain using the mean-field (A^ = oo) result for the 
correction term 



Xu{T) 



0.86X + 0.145- 



knT 



(5.3) 



C. Dynamic susceptibility and structure factor 

In this subsection, we compute the scaling functions ^is{k,uj,x = oo, yp = 0,y^ = 
0) = ^s{k,'oJ) and Ei{k,ZJ,x = oo, yp = 0,y^ = 0) = E{k,Lj) for staggered dynamical 
susceptibility and structure factor at the critical point g = Qc- These two scaling functions 
were introduced in ( [1.26| ) and are related by the fluctuation-dissipation theorem S(A;,ZI7) = 
lm.^s{k,uj). As before, ^s{k,uj) is simply related to the polarization operator for z— fields: 
$s(/c,a;) = {kBT/2{hc±)'^) Il{k,uj). The limiting behavior of 11 and hence $^(^5^) at small 
and large k and uJ can be obtained by properly expanding the real part of ( |4.13|) . For large 
k, UJ, we found using the results of P5| , |5B| 

^sik,uj) = , -- + ^^ — + 5.4 



On the other hand, at small momentum and frequency, an expansion in the real part of 
( TO ) yields 



VE I A;^(l + 26/75) 



Re ^s{k,u) = - '^^^'720^^ " + O [{k,uY) 1 (5.5) 

We consider next 1/A^ corrections to these results. At small k and ID, the expansion in 
1/A^ does not involve logarithms. Regular 1/A^ corrections to $s were found to be quite 
small for unfrustrated antiferromagnets ||2^ and we expect the same to be true in our case 
as well. On the other hand, at k,uJ ^ 1, the behavior is nearly the same as at the critical 
point at T = 0, and using the results of Appendix C we found that the leading term in ( p.4|) 
is modified to 

*^(^'^)=16(g^-(J^+z5)2)i-'^/2 ^'-'^ 



where Aj^ = 1 + 0{1/N) and fj is given by 

It is also not difficult to compute explicitly the imaginary part of the polarization op- 
erator, which then yields the scaling function for the dynamic structure factor. In the two 
assymptotic limits of large and small ZJ we obtained 

A.sin(.,/2) ejUJ^-k') 
^ ' 16 ^2)1-^/2 ^ ' 
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for io ^ 1 and 

for tJ ^ 1 and /c ^ 1. In (^.7|), 6{x) is a step function. It is also not difficult to obtain 
the N = oo expression for S for both oJ <^ 1 and /c -C 1, but in this region of momentum 
and frequency, quasiparticle excitations are overdamped, and one again cannot restrict to 
the N = oo result of collisionless Landau damping. We can only expect on general grounds 
that at small uJ and arbitrary k, a;) oc uJ. 

We consider, further, the static structure factor S{k) defined by ( |4.19| ). In the quantum- 
critical region, we have at g Qc 




m = ^ I - I I I Hk) (5.9) 

where I{k) = J duo (1 — e^"^)^^ 'E{k,uj)/'K. Notice that the functional form of S{k) is 
similar to that for unfrustrated antiferromagnets ||25|. Moreover, in both frustrated and 
unfrustrated cases, E{k,ZJ,oo) behaves at large frequences as E{k,uj, oo) oc (cZ;)"^^''. The 
difference between the two cases is in the value of f/. For unfrustrated antiferromagnets, 
f] ^ and the frequency integral in I{k) is convergent. For frustrated systems, f/ > 1 at 
least, at large (see ( p.4|) ), and the integral over frequency in I{k) is divergent, which 
actually means that the dominant contribution to S{k) at small temperatures comes from 
the frequencies of the order of a cutoff. Specifically, using ( |4.13| ) we obtain 



Uk) = [fA^-i - 6^-^)1 + I'(k) (5.10) 

?7 — 1 Lv n 

Here = 1 + 0{1/N), 6 was defined after (|5.1|) , A = hcA/kBT where A is a relativistic 
cutoff in the theory, and I'{k) is a universal function of momentum, which as /c tends to 
/'(O) = 1.67 + 0(1/A^). The nonuniversality in I{k) at low temperatures is now transparent. 
Recall however that ( |5.1U| ) is valid only in the quantum-critical regime where the temperature 
is the only scale for fluctuations. In the renormalized-classical regime, the analog of I'{k) is 
proportional to the square of the actual correlation length ^, and is exponentially large at 
low T compared to the nonuniversal piece in I{k), which does not contain any dependence 
on ^. 



D. Local susceptibility and spin- lattice relaxation 



Unlike S{k), the momentum- integrated dynamic structure factor Sl{uj) = 
J (fk S{k, uj)/A'k'^ is universal in 2d, as we already found in the renormalized-classical region. 
The scahng function for Siiyj) was introduced in ( p..36| ). At g = gc and N = oo, this scaling 
function can be deduced directly from ( [4.13| ). A simple calculation yields 



Kip) 



327r 



21og- 



1-e-® 



+ e(cJ-20) cJ-2e + 21o^ 



l-e-0 



(5.11) 
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where K^u) = Ki{u, xi = oo, ,yp = 0, = 0). At small u, this reduces to 



while at ^ 1, we have 



KP) - ^ (5.13) 



Note that linear dependence on uJ is present in both limits (in the unfrustrated case, K{uj) 
saturated at large a;). 

We now consider corrections to these results. At cJ ^ 1, the expansion in is free 
from divergences because each z— field propagator has a gap O ~ 1. The expansion in 
then holds in integer powers of and numerically we expect the corrections to ( [5.12D to 
be quite small. On the contrary, at large cZ; 1, the actual form of K(ZJ) is different from 
the N = oo result because of the singular corrections. Using ( |5.7| ) we obtain instead of 

^ ' 327C f] ^ ' 

Finally, the cJ — ^ limit of Sl{uj) is related to the transverse spin- lattice relaxation rate. 
As before, we assume that the hyperfine coupling constant tends to a finite value at = 
and the dominant contribution to 1/Ti thus comes from the momentum range k = 0(1). 
Using ([OGp , i ^M ) and (CT) , we then obtain 



(5.15) 



where Z = {Vb - 1)/8N + 0(1/A^2)^ 

A parallel analysis can be done for longitudinal spin-lattice relaxation I/T2 defined in 
( [4.28| ), and the result is 




Nl_h fNksTy-' 



(5.16) 



VI. QUANTUM DISORDERED REGION 

Let us first notice some crucial properties of the quantum-disordered {g > Qc) phase at 
T = 0. The presence of free spin-1/2 z quanta implies that Xs{k,uj) only has a branch cut 
in the complex uj plane. This should be contrasted with the behavior of coUinear antiferro- 
magnets |2^ which had an additional spin-1 quasiparticle pole. We will compute below the 
structure of this branch cut at iV = 00. 

Before doing this, it is useful to introduce our precise definition of the prefactor A. 
We will use the T = form of the local dynamic structure factor Sl{ijj) to specify the 
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normalization. Using the fact that the dynamic susceptibihty involves a response function 
with 2 spinon intermediate states, each with a gap A, and that each spinon propagator has 
a quasiparticle pole we can show quite generally (to all orders in 1/N or e = 4 — D) that 
near the threshold in the quantum disordered phase we must have 

Sl{uj) = A ^^ ~^^ e{huj - 2A) u close to 2A. (6.1) 

The above form defines the values of A and A. Combined with c, these parameters univer- 
sally determine the entire staggered susceptibility. 

The N = oo computation of Xs is standard ||25|. The z propagators acquire a gap A 
given by 



A = 47r (6.2) 




We evaluated the susceptibility using ( |1.11|) and found at = oo 

Imx.(A:,^) = sgn(.;)::^-^_==^^ {u' - c^k' - AA^h') (6.3) 



with A = g'^ A / {AT{h^) . These results are consistent with (|1.28|) and ( |1.32| ) provided z/ = 1 
and f] = 1. Note that Xsik,^) has branch cuts emanating from ±(4A^//i^ + c^k'^) to ±cxo. 



Compare this result with the confined spinon model of ( |2^) where at = oo, Im^^ was 
simply a delta function. 

It is simple to extend the above results for Xs to finite temperature and to order 
The main effect of small T is to fill in the gap in the spectrum by exponentially small 
terms. The corrections do not introduce any essentially new features, and will not be 
considered here. 



VII. APPLICATION TO A S = 1/2 ANTIFERROMAGNET 

In this section, we compare our scaling results to the properties of the S = 1/2 Heisen- 
berg antiferromagnet on a triangular lattice. As input, we need the values of sublattice 
magnetization, two spin stiffnesses and two spin-susceptibilities at T = 0. In Appendix B, 
we have calculated these parameters in the 1/5* expansion, to order 1/S for the stiffnesses 
and susceptibilities and to order 1/5*^ for sublattice magnetization. The extension of our 
large S results to 5* = 1/2 yields 

No = 0.266; x± = ^ X = Ps = 0.086J, c = = l.OlJa (7.1) 

The magnitude of the 1/5*^ result for the sublattice magnetization indicates that higher-order 
corrections are rather small. 

Let us now summarize the scaling predictions which follow from the values in ( |7.1| ) . For 
the uniform susceptibility, we obtain 
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\ ha J J 



0.084 + 0.07- 



J 



(7.2) 



in the renormalized-classical regime, and (using the N = oo result for the correction term) 



Xu 



\ ha J J 



0.072 + 0.14- 



J 



(7.3) 



in the quantum-critical regime. Comparing ( [7.2|) and ( [7.3|) , we observe that the slope of 
Xu in the quantum-critical regime is nearly twice as large as in (|7.2|) , while the value of 
the intercept is larger in the renormalized-classical regime. Further, the correlation length 
behaves in the renormalized-classical regime as 



e~0.24(^-^j exp[A7TpJkBT] 

where 47rps ^ 1.08 J, and deep in the quantum-critical region as 

^ O.SlJa 



(7.4) 



(7.5) 



Finally, in the renormalized-classical regime, the universal contribution to S{k) is dominant, 
and for k = we obtain from ( [4.16| ) and ([4.19|) 



(7.6) 



In the quantum-critical region, the dominant piece in 5(0) is a temperature-independent 
contribution from lattice scales, and we can only conclude that deep inside quantum-critical 
region, S{0) ~ y4 + B{T/Tq)^~^, where A is a T -independent nonuniversal piece. Using the 
large N results we found B ^ — 0.27a^ and Tq ^ 0.54J. 

From the discussion in the bulk of the paper, we expect the crossover between the classical 
and quantum regimes to occur somewhere around xi = 1 i.e. at ksT = 27ips ~ 0.5J. This 
indeed is not a very small crossover temperature. However, the analysis for the unfrustrated 



case 



^ shows that the uniform susceptibility displays quantum-critical behavior starting 
already below Xi = 1. We therefore first compare our results with the numerical data on 
uniform susceptibility. 

The temperature dependence of Xu was recently studied in high-temperature series ex- 
pansions for 5* = 1/2 triangular antiferromagnet The data show that Xu obeys a Curie- 
Weiss law at high T, passes through a maximum aX T ^ 0.4J, and then falls down. In gen- 
eral, the temperature where Xu has a maximum roughly separates the low-temperature region 
below the maximum where a long-wavelength approach is valid, from the high-temperature 
region where the physics is dominated by lattice-scale effects. It is unfortunate that this this 
temperature is rather low for the triangular antiferromagnet, because it reduces substan- 
tially the temperature range for low-energy behavior (for comparison, in the square-lattice 
antiferromagnet, the maximum in Tc occurs at fcsT ~ J). Numerical data is available 
only over a small T region below the maximum. Nevertheless, we fitted the data by a linear 
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in T dependence and found 0.13 ±0.03 for the slope and around 0.06 for the intercept - both 
results are in reasonable agreement with our quantum-critical expression (fTSl) . Finally, at 
very low T, we expect a crossover to the renormalized-classical regime, and the T = value 
in (|7.2|) is also consistent with the data. 

Now about 5'(0). Previous studies of square-lattice antiferromagnets have shown that 
we can hardly expect to observe pure quantum-critical behavior for 5(0) at xi ~ 1. Indeed, 
the leading correction to 5(0) due to the deviation from purely quantum-critical behavior is 
55(0) = C{27ips/kBT), where at = cxd we found C = 0.45a^. Clearly then, at ksT ~ 0.5 J, 
temperature dependence related to deviations from pure criticality is likely to overshadow 
the weak ((T/Tq)''"^) temperature dependence in 5(0) at g = g^, instead, we expect that 
at such T, the structure factor should roughly follow 5(0) = A + C{2TTps/kBT), or (still 
considering second term as a correction) ksT log S{0) ~ 2ttPsC/A + ksT log A. The series 
expansions |^0| yielded ksT log S{0) which increases linearly with T upto about 0.5 J. This 
is consistent with our crossover expression, but inconsistent with the renormalized-classical 
formula, ( |7.6| ) which predicts that ksT log S{0) decreases with temperature. We therefore 
do not believe that the numerical data correspond to the classical regime, as was suggested 
in Ref |0|. 

Finally, the correlation length. Series expansions reported that ^ is approximately one 
lattice spacing at fc^T = 0.4J. This is substantially lower than our renormalized-classical 
result ~ 5a at the same temperature, but is consistent with the value of ^ deep in the 
quantum-critical regime C,quant ~ 1.25a. We emphasize however that Ref. [0 defined ^ as 



= —{1/ S{k)){dS/dk'^)\k=o - this definition yields a nonuniversal value of ^ for quantum- 
critical frustrated antiferromagnets. On the contrary, our definition of ^, from the long- 
distance decay of spin-spin correlator, always yields a universal result. Besides, even if 
the universal piece in S{k) is dominant, as in the renormalized classical regime, the two 
definitions are still nonequivalent even at = oo simply because spin structure factor 
is related to the polarization operator of z— fields, which unlike field propagator, does 
not have a Lorentzian form. In the classical regime, the rescaling factor between the two 
definitions of C, is ^series = (2/3)^ours at = cxd. The value of the rescaling factor in the 
intermediate and quantum-critical regimes is difficult to estimate, but on general grounds 
it should be smaller than 2/3 because the nonuniversal piece in 5(0) becomes dominant at 
g = gc- We therefore expect that the actual correlation length is in fact larger than reported 
in This again is consistent with our observation that at /c^T around 0.4 J, the system 
is in the crossover region between renormalized-classical and quantum-critical regimes, and 
is probably closer at ksT = OA J to the quantum-critical regime. 



VIII. CONCLUSIONS 

In conclusion, we summarize our main results. We have presented a general scaling frame- 
work to describe frustrated antiferromagnetic systems near the quantum phase transition 
between classically ordered and quantum-disordered ground states. 

We considered various scaling functions for experimentally measurable quantities both 
on the ordered and disordered sides of the quantum transition and have shown that the 
observables which probe the behavior of antiferromagnets at low energies are completely 
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universal functions of just a few measurable parameters at T = 0. On the ordered side, 
these parameters are sublattice magnetization and transverse and longitudinal spin- stiffness 
and spin susceptibility. 

We then specialized to particular field-theoretic model of the transition (results for a 
different model are briefly noted in Appendix A) Our approach began with the fundamental 
assumption that the disordering transition at T = is continuous and that vortex-like 
excitations with a nonzero local Zi flux are irrelevant at low energies. We showed that, in 
this situation, the proper low-energy theory near the transition is given by the SU{2) x U{1) 
sigma-model for spinon fields. All physically observable excitations are collective modes 
of two spinous. The global SU{2) symmetry of the sigma-model action is related to spin 
rotations, while the global U{1) symmetry is related to lattice transformations |2^. For 
triangular and other commensurate noncoUinear antiferromagnets, this lattice symmetry in 
fact reduces to a discrete symmetry (Z3 symmetry in case of triangular antiferromagnets). 
At the quantum transition point, the symmetry of the action enlarges to 0(4). 

We then extended our action to a general by considering spinous as A^-component 
objects, and used the powerful technique of 1/A^ expansion. The extended action has 
SU{N) X f/(l) symmetry. The fixed point in this approach has its internal symmetry 
enlarged from SU{N) x f/(l) to 0{2N) for any N. 

We then used the 1/A^ expansion to explicitly compute the scaling properties of the field- 
theory, always finding that they were consistent with the more general scaling ansatzes. We 
made definite predictions for the dynamic structure factor, static susceptibility, correla- 
tion length, local and static structure factors, and the spin-lattice relaxation rate in the 
renormalized-classical and quantum-critical regions. We also briefly discussed the low-T 
behavior in the quantum-disordered region. 

Finally, we compared our results to the properties S = 1/2 triangular antiferromagnets. 
We determined the input parameters in the scaling function from a. 1/S expansion on the 
original lattice Hamiltonian, and made quantitative predictions about the form of uniform 
susceptibility, correlation length and static structure factor. We compared the results with 
the data of recent high-temperature series expansions. All of the data were consistent 
with the interpretation that there is a narrow window of quantum-critical behavior just 
below the temperature at which the uniform susceptibility passes through its maximum. 
However, more detailed numerical and experimental results are needed before any definitive 
conclusions can be reached. We hope that it will be possible to perform measurements 
in a T range between the 3d ordering temperature and the temperature where uniform 
susceptibility has a maximum. Our prediction is that in between the two temperatures, the 
uniform susceptibility should follow our formula for the quantum-critical regime. 
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APPENDIX A: FIELD THEORY WITH CONFINED SPINONS 



In the event there is continuous transition from the magnetically ordered state to a 
quantum-disordered state with confined spinous, we expect that it can be described by a 
continuum field theory of the rii and 112 fields themselves. All fields are now singlets under 
the Z2 gauge symmetry and the Z2 vortices are permitted. Such large-M field theories have 



been considered earlier by Kawamura and Azaria et. al. |29]. A potential problem with 
this approach is that the results of the D = 2 + e analysis |T^ are not obviously consistent 
with the D = A — e and large M theories. The universal properties of such nearly-critical 
antiferromagnets are rather similar to those of the coUinear antiferromagnets considered in 



Ref. ||25|. Therefore we will be rather brief, as the analog of all the results in the body of 
this paper can be obtained by minor modifications of those of Ref P5 . 
We will consider the action 

(Al) 

where pi^^ = p*}^^, P2,^t = (P||^ — 2p° ^^)/4. The potential V^(ni,n2) can either impose the 
hard-spin constraints (in a D = 2 + e expansion) 

= = 1 ; ni ■ n2 = (A2) 

or the soft-spin potential (in a = 4 — e expansion) 

V{ni, na) = ^ro (n? + ^4) + Ui (n? + + U2 (ni x 112^ (A3) 



Kawamura [^] introduced a large M expansion of ([Al|) in which the vectors rii, 112 are gen- 
eralized to M-components; the action then has a 0{M) x 0(2) symmetry. The relationships 
betwen the large M, e = D — 2, and e = 4 — D expansions have been discussed by Azaria 
et. al. PI. 



Here we will discuss some simple properties of the large M expansion. The results 
have striking differences from the large N expansion of this paper, in particular, the phase 
transition at M > 3 belongs to the universality class different from 0{M + 1) model. Which 
of these two expansions is more appropriate for the physical case M = 3, = 2 is not 
quite clear, and numerical studies of frustrated antiferromagnets will be quite useful in this 
regard. The most obvious difference is of course in the absence of spinous in the large M 
theory. The staggered susceptibihty Xs{k,u!) now has delta-function quasiparticle peaks, in 
contrast to the branch cuts of the large theory. Differences also appear in the behavior of 
the correlators of the conserved charges and currents. A key property of the 1/M expansion 
is that the p2,/x couplings are irrelevant. This immediately implies that the universal ratios 
of the stiffnesses (Eqn (|1.21| )) obey Tp = = 2 at M = cx). We computed T to first order 
in 1/M and obtained 



2fi 

T„ = = 2 ^ (A4) 
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In performing the 1/M calculations, we introduced three Lagrange multipliers in the func- 
tional integral to impose the constraints (|A2|) , and also introduced condensates of the rii 
and n2 fields. The computations are a bit tricky: we found that the 1/M correction to T 
is related to the difference in the Green functions of the transverse components of ni 2 and 
the fluctuating components along the directions of the condensates. This difference clearly 
disappears ai g = Qc] however the correction to the ratio of the stiffnesses (i.e., to to O(fc^) 
in the full Green functions) remains finite at g = Qc because it includes integrals which are 
divergent aX g ^ Qc- 

Also interesting is the behavior of the uniform susceptibility Xu{T) in the quantum- 
critical region. It is simple to show that at M = 00 this is given by precisely twice the 
mean-field result of Ref [^], with c ^ c±. Contrast this with the result of the deconfined 
spinon model of the body of the paper: there we found that the N = 00 result was one-half 
the result of Ref [EH] ! 



APPENDIX B: SPIN- WAVE CALCULATIONS AT T = 

For experimental comparisons of the results obtained within 1/N expansion, we need 
the T = expressions for sublattice magnetization, spin-wave velocities and uniform spin 
susceptibilities. Below we will calculate these quantities for the Heisenberg antiferromagnet 
on a triangular lattice in an expansion in 1/25, where 5* is the value of the spin. Though 
we will use large S approach, our chief interest is in the case of S* = 1/2 when quantum 
fluctuations are the strongest. As we will see below, the convergence of the perturbative series 
in 1/25* in triangular antiferromagnets is very good (as it is on the square lattice 
and the 1/5 expansion is likely to give quite accurate values of observables, even for 5 = 1/2. 

We now turn to a description of the calculations. We consider here the model with 
interactions between nearest neighbors: 

7^ = jJ^SiSi+A. (Bl) 

The procedure of doing the 1/S expansion is rather standard and involves several steps 
which include (i) the transformation from spin operators to bosons via Holstein-Primakoff, 
Dyson-Maleev, or some other transformation, (ii) the diagonalization of the quadratic form 
in bosons, and (iii) the use of a standard perturbative technique for Bose-liquids to treat 
the interaction between spin waves. Noninteracting spin waves have energy which scales 
as S, while the interaction vertex involving m bosons scales as S'^"''"/^; this gives rise to 
an expansion in powers oi 1/S for anharmonic contributions, similar to that in a weakly 
interacting Bose gas. 

Another important issue related to the 1/S expansion, is the number of bose fields which 
one has to introduce in order to keep track of the whole spin-wave spectrum, not just the 
low-energy modes. This is important because quantum fluctuations are not divergent in 2d, 
and the 1/S expansion involves sums over the whole Brillouin zone. In the general case, the 
number of different bose fields is equivalent to the number of magnetic sublattices. However, 
in several special multisublattice magnetic configuration can be transformed into a 

one-sublattice ferromagnetic one by applying a uniform twist on the coordinate frame. In this 



35 



situation, the spin- wave spectrum has no gaps at the boundaries of the reduced Brillouin zone 
and one can describe all excitations by a single bosonic field, as in the case of a ferromagnet. 
Obviously, the triangular antiferromagnet in a zero magnetic field is an example of such 
special behavior: the 120° ordering becomes a ferromagnetic one in the twisted coordinate 
frame with a pitch Q = (47r/3, 47r/v^). We therefore will use a one-sublattice description 
of triangular antiferromagnet whenever possible. This indeed substantially simplifies the 
calculations. 

We start with the transformation from spin operators to bosons. The choice of the 
transformation is indeed only a matter of convenience, and the final results are indepen- 
dent of the way how bosons are introduced. Nevertheless, there are several possibilities 



extensively discussed in the literature We found it most convenient to use here the con- 
ventional Holstein-Primakoff transformation because it preserves the Hermitian properties 
of the Hamiltonian. We therefore use 



Sz = S-a^a- S^^ = \ 2S - a^a a; S' = a\ 2S - a^a (B2 



Substituting this transformation into ( |B 1| ) , expanding the radical, and restricting to only 
cubic and quartic anharmonic terms, we obtain after some algebra 

n = no + 3JS{n2 + n3 + n^) (bs) 

where Ho = ~\JS'^N is the classical ground state energy, and other terms are 

k ^ 

"^4 = H ai4«3a4 [4(z/i_3 + v^-^) + z/i + z/g + z/3 + z/J 



-2 (^a| 020304 -I- 04030201) (z/i + z/2 + z/3) 



Here i = ki, and 



/ 3 

^3 = M/ -377 yi(44«3 - ala2ai){i?i + U2). (B4) 
V 00 



1 / J k^ ky^/^^L _ 2 , kx ( kx ky\^^\ /T-ir\ 

i^fc = 2 I cosfc^ + 2cos-^cos^^ I ; z/^ = - sm I cos y - cos I . (B5) 

Finally, Ak and Bk are given by 

ly 3 

^ = l + y; Bk = --Uk (B6) 

At S* = 00, anharmonic terms are absent and Tii can be diagonalized by a standard Bogol- 
ubov transformation 

ak = lkick + Xkclk) (B7) 
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with 



Xk = -r-^ -r ■ (B8) 



2Ek J \Bk\ \Ak + Ek 

and 

Ek = (Al - Elf I' = ((1 - Uk)il + 2vu)f" (B9) 
The diagonahzation yields 

ni = Y.Ekcick (BIO) 



It follows from Eqn. ( |B9| ) that the excitation spectrum of the ideal gas of magnons has 
three zero modes, as it indeed should. Two of these modes are at = where Q = 
(47r/3, 47r/v^) is the ordering momentum in triangular antiferromagnet, and the third is at 
k = and describes soft fluctuations of total magnetization. The expansion near zero modes 
gives two spin-wave velocities 

c± = c±Q = ^-^JSa (Bll) 



Cfc=0 



3^3 



JSa 



(B12) 



The ratio of the two at 5* = oo is cy /c± = V2. This was also obtained in other ap- 



proaches 1 26 



The infinite 5* spin-wave results can be also used to get the first quantum correction 
to on-site magnetization [^. Indeed, (a^a) in ( [B2D is nothing but the density of particles 
which is finite due to the anomalous term in the quadratic form. From (|B7| , p8| ), we have 
{a\ak) = {Ak — Ek)/2Ek, and therefore noninteracting spin waves reduce the sublattice 
magnetization to 



<S>=S 1 



2S 



E 



Ak — Ek 
Ek 



SI 



0.522 
2S 



(B13) 



We next consider corrections to Eqns ( piO|) and (pi2|) due to the interactions between 
spin-waves. We will follow the same line of reasoning as for square-lattice antiferromagnets. 
However, the presence of cubic terms makes the analysis considerably more involved. 

We start with the spin-wave velocity renormalization. 



1. Spin- wave velocity 

Our goal is to obtain the leading 1/S renormalization of spin- wave excitations. For this 
we consider first-order self-energy corrections due to quartic anharmonicities and second- 
order corrections due to cubic anharmonicities (recall that cubic terms have the overall 
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factor S*^/^). The corrections due to quartic terms are easy to compute, because to leading 
order in 1/5*, one can get away with simple one- loop diagrams. Equivalently, one can simply 
decouple the four-fold term in eq. ( p^) by making all possible pair averaging. The quadratic 
form allows for nonzero normal {a\ak) and anomalous (aka-k) pair products of Bose particles, 
and the decoupling changes At and Bk to 




2 2S 2S^ Ep \ A "Jj 8S^ Ep ^ ' 

A simple inspection then shows that the renormalized spectrum [E^ = {Al — still 
keeps a zero mode at k = 0, but acquires a finite gap at k = ±Q: 

^2^_AyMLzi^ (B16) 

8S ^ Ep ^ ' 

This finite gap is indeed an artifact of using only quartic terms, and cubic anharmonicities 
should restore the correct structure of the spectrum, as we demonstrate below. 

There are several ways to deal with the cubic terms: one can either calculate the effective 
four- fold vertex produced by two triple vertices p5| , |iD[] , and then use the decoupling proce- 



dure, or one can transform to quasiparticles (i.e., diagonalize the quadratic form) considering 
first only quartic corrections, and then calculate the renormalization of the excitation spec- 
trum due to cubic terms in the second-order perturbation theory. Below we use the second 
approach which is technically advantageous. We therefore first transform from particle op- 
erators (ofc) to quasiparticles (c^) using eq. (|B7|), but with Ak and Bk instead of Ak and 
Bk- The bare Hamiltonian then keeps the form of eq.( |B10|) with Ej^ instead of Ej^. On the 
other hand, the structure of cubic vertices becomes more involved after the transformation 
to quasiparticles, and instead of Eqn. (p^) we obtain 



^'"iJs ^ cI4c3<fi(l,2;3) + yAclMl,2,3) + H.c (B17) 
The vertex functions $i and $2 are given by 

*.(l,2;3) = *liiM; *,(1,2,3)=*#M (B18) 



where 



E1E2E3 J E1E2E3 



$2(1,2,3) = z>,/«(/f /f - /l^Vi^^) + ^2/l^^(/f /f - /l^Vi^^) + - /i^Vi^^; 

and 
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The self-energy diagrams to order 1/5* are shown on Fig. ^. We see that cubic terms give 
rise to both normal and anomalous self-energy parts so that the dispersion relation again 
has the form typical for a 2 x 2 problem: 

{u + uj)f = (Ek + J^sik, iu)f - (B20) 

where T,s^a{k,uj) = | iT,^_{k,Lj) ± T,^_{—k, —u). However, it is not difficult to check that 
S_ _ ~ ^+,+ ~ 1/5* and therefore anomalous self-energy terms contribute to the excitation 
energy only to order l/S*^, while to order 1/S a. solution of Eqn ( P20|) is simply uj = 
where 

El = El + 2EkJ:+,.ikEk) (B21) 

We therefore need to evaluate here only the normal component of the self-energy. The 
analytical expression for _ is 

^ 3 / |$i(l,2;fc)p |$2(l,2,fc)n 

""-"'-^^^ ""'^ --^S^ Ui + E,-E,^ E, + E, + E,] ^^22) 

To leading order in l/S" we can indeed use nonrenormalized values for Ak,Bk, E^ in the r.h.s. 
of ( p22|) . 

We first demonstrate that E^ has a true zero mode dX k = Q. For this we need to 
evaluate E+ _(Q,£'q). We found the following equahty to be quite useful in the calculation 

^q±Q/2 = {Aq±Q/2 + Bg±Q/2) - - (B23) 

Substituting ( P23|) into the expressions for the vertex functions and using Aq = Bq = 3/4, 
we obtain after simple algebra 

$,(l,2;Q) = $,(l,2,Q) = i^^l^ (/^Vf -/i^Vi^^) (B24) 

Substituting, then, the vertex functions into the formula for the self-energy we obtain using 
( [BT9D 

Finally, upon substituting this result into Eqn ( [B21| ) and using ( [B16| ) for Eq, we find that 
the gap in the excitation spectrum disappears as it should [pO[] . 

Our next step is to expand E^ and near the zero modes, and obtain the corrections 
to the spin- wave velocities to order 1/5*. The expansion near = is quite straightforward 
because $i(l,2;/c) and $2(152, /c) both scale as k at small fc, and one can therefore safely 
neglect E^ in the denominators in ( |B22| ). Doing the algebra, we obtain the renormalized 
spin-wave velocity at /c ~ in the form 
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1 

2S 



1 

35 



Eh 



Q 



(9 - 4QD 



(B26) 



where 



sm h sm + sm — 

2 4 4 



1 
3 



sm 



4 "^a; 



sm + sm 



IQl 



Numerical intergation then gives 



0.115 
2S 



(B27) 



The structure of the expansion near k = ±Q is more involved and we refrain from presenting 
the analytical expression for the spin-wave velocity. Numerically, we obtained 

Comparing ( [B27| ) and ( [B28D , we observe that quantum fluctuations tend to diminish the 
difference between the two spin-wave velocities. This is consistent with our result ( p.7|) that 
the relative difference between c± and cy should disappear at the quantum-critical point. 
We will use ( [B27| ) and ( |B28| ) below and now proceed with the calculations of sublattice 
magnetization. 



2. Sublattice magnetization 

We have shown above that to leading order in 1/S, the correction to sublattice magne- 
tization comes already from noninteracting magnons (Eqn ( |B13|) ). Here we obtain the next 
term in the expansion in 1/S, which is also the leading 1/S correction to the density of 
particles. We again have to consider both quartic and cubic terms, since they contribute at 
the same order to J2k{(4:^k)- As before, quartic terms only renormalize the coefficients in 
the quadratic form, and hence change the expression for the density of particles to 

k I I hk 

where Ak and Ek are given by ( |B14| ) and ( |B15| ). In explicit form 

1 ^ ^ _ 1^ 1+ Z/fc/2 _ _9_ ^ ^ _ ^ Z/p(l - Up) ^ yq{l - Uq) 

2 ^ -Efc 2 ^ Ek 325 p Ep q Eq 32S p Ep ^ E^ 

(B30) 

We see that the very last term behaves near q = Q as \q — Q\^^ which makes the integral 
over q divergent. The divergence is indeed an artificial one and should disappear when we 
add the contributions of the cubic terms. 
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To see how cubic terms modify ( [B29| ), we express the density of particles in terms of the 
quasiparticles operators using (p7| ) and (p^ ): 



E(4«.) = -^ + ^e|^-e|^ <ckc^'^> + <cW> (B31) 

The first two terms are just the renormahzed spin- wave terms. The third correction is 
related to the anomalous self-energy term in Fig.^. Performing the frequency summation in 
this term, we obtain 

where 

^'-^E^E, E, + E, + E, (^22) 

Finally, the last term in (|B31| ) contains the density of quasiparticles. This density is finite 
to order 1/ S because among cubic non-linearities, there is the term which describes simulta- 
neous emission of three spin-waves. Evaluating the expectation value of (c|.Cfc) by the usual 
means, we obtain 

, El ^ '''' ? ' ^ ^ 

where 

T, = J:J— (B35) 
E1E2 {El + E2 + Ek)'^ 

We first show that the total expression for the density of particles is free from divergencies. 
Simple inspection of Eqns ( |B32| ) - ( [B35| ) shows that the divergent contributions from the 
cubic terms (namely, 1/E'^ and l/E"^ terms in ( P32|) and 1/E'^ terms in ( P34|) ) come from 
the region k ^ Q, where \l/ and T tend to constant values. For these k, we again use ( |B23D , 
substitute it into the vertex functions, and after a simple algebra obtain 

= -(1 - ^q) E ^^^^^^ - E,Tq + 0{El) (B36) 
P Ep 



Substituting further this expression into ( |B32D and comparing the result with the divergent 



piece in ( P30|) , we find that the 1/E^ contributions from cubic and quartic terms, and the 
l/E"^ contributions from the two cubic terms cancel each other, so that the 1/5* correction 
to the density of particles is finite, as it of course should be. We then performed numerical 
computation of the 1/S terms in (|B31|) and obtained 



/ n n27\ 

(S) =[S- 0.261 + (B37) 
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For S = 1/2, Eqn ( P37|) yields (S) ~ 0.266, which is close to half of the classical value. A very 



similar result was obtained earlier by Miyake |51[], who calculated the on-site magnetization 
to order 1/5"^ by evaluating numerically the response to a staggered magnetic field. His 
estimate for the l/S"^ correction is however somewhat smaller than ours (0.01 instead of 
0.027). In any event, l/S*^ terms are rather small and can hardly change substantially the 
lowest-order spin- wave result for the magnetization We therefore found no support for 
the recent claim |^ that the value of magnetization is substantially lower than the spin-wave 
prediction. Note, in passing, that for square lattice antiferromagnet, the first anharmonic 
correction to {S) is identically zero |^ . Indeed, cubic terms are absent in the square-lattice 
antiferromagnet, and 1/S corrections due to quartic terms do not change the shape of the 
quasiparticle spectrum (that is, Ak/Ek = Ak/Ek). The next to leading order correction 



in the square-lattice case has been calculated and found to be very small |52 



3. Uniform susceptibility 

Now we calculate, to order the response of a triangular antiferromagnet to an 

external magnetic field. We have already discussed in Sec.|ITI|that the magnetic susceptibility 
tensor in a triangular antiferromagnet has the form pi 



Xa p = X^^ap + (X|| - Xi.)mamp (B38) 

where m is a unit vector which specifies the plane of spin ordering. This form of Xap implies 
that the ordered state should have two different spin susceptibilities. They can be viewed 
as the response to the field applied perpendicular to the plane of spin ordering, i.e., along 
m (x||)) and as the response to a field directed in the plane (x±)- In the latter case, we need 
to introduce an infinitesimally small anisotropy which keeps the spins in the basal plane. 

For classical spins, the transverse and longitudinal susceptibilities can easily be obtained 
by minimizing the ground state energy. This yields x± = X\\ = 2/Q\/?)Ja? where a is the 
interatomic spacing (a^\/3/2 is the unit cell volume). As in the bulk of the paper, we define 
Xi_ and x\\ without the gyromagnetic ratio gfis/fi'- We see that the two susceptibihties are 
equal in the classical limit . This degeneracy in the response to a magnetic field in a 2D 



triangular antiferromagnet has attracted some attention in the past as an example of the 



'order from disorder" phenomenon [^,^,^,^,0. For our present purposes, it is sufficient 



to observe that the degeneracy is a purely classical effect. It is not related to the symmetry 
properties of a quantum system and therefore should be broken by quantum fluctuations. 

Technically, the computations in a finite field are more involved because the transverse 
field breaks the 120° ordering in the basal plane. In this case, a transformation to a twisted 
coordinate frame is no longer advantageous because umklapp processes also contribute to 
order 1/S. It is then more convenient to introduce a separate bose field for each of three 
sublattices. For the longitudinal response, the 120° ordering in the basal plane is preserved 
and a one-sublattice description with no umklapp terms is still valid. However, one has to 
be careful in this case as well, because in the presence of a field, the excitation spectrum is 
no longer an even function of k. This is consistent with the fact that time reversal symmetry 
in a magnetic field requires that in changing k — >■ —k in the spectrum, one has to change 
simultaneously the sign of H. 
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The corrections to the susceptibihty tensor to order 1/ S were computed by Golosov and 
We refrain from presenting the details of the calculations and list here only 



one of us [56 



the results. To order l/S", they are (notice that the definitions of x± and x\\ iii ^ire 
interchanged compared to ours): 



where 



and 



1 

2^ 



9V3Jc 



1 

2S 



E 



Eh 



E 

k 



3 

2^ 



XII 



9V3Jc 



Eh 



E 



0.448 
2S ■ 



(1) /(2)^ 



f'-' f- 

Eki + -E'fca -E^fci Ek2 



(B39) 



(B40) 



0.291 
2^ 



(B41) 



where /'•*•' = f^''^'^ were defined in ( |B19D . Note that contrary to the situation in a stacked 
3d triangular antiferromagnet where x\\ > X±) the transverse (in-plane) susceptibility in 
the 2d case turns out to be larger than the longitudinal one; this gives rise to an uncon- 
ventional phase diagram in a magnetic field which has been discussed several times in the 
literature 



4. Spin stiffness 

With the values of the two spin-wave velocities and spin susceptibilities at hand, we are 
now in a position to calculate the spin stiffnesses. To order 1/S they are 

P± = X±cl = ^JS'Zi pii = xiicj = ^JS'Zl (B42) 

where 

= 1 - 0.125/25, = 1 - 0.678/25 (B43) 
Finally, substituting the results ([B39D- ([B42[) into (|2.6|), we ontain for = 2 



APPENDIX C: COMPUTATIONS IN THE NEEL STATE AT T = 

In this Appendix, we derive to order 1/A^, the expression for T = sublattice magne- 
tization near the quantum phase transition. This result will be used in the derivation of 
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the universal scaling forms for uniform and staggered susceptibilities in both renormalized- 
classical and quantum-critical regions. We also reproduce the expressions for the relative 
difference between longitudinal and transverse spin-stiffness and spin susceptibility which 



were obtained by other means in Sec |ITl 



Our point of departure is the functional integral for the SU{N) x U{1) a— model, Eqn 
At T = 0, the spin-rotation symmetry is broken, and we represent the N component 
complex vector z of length as 

z = {ao + iai, tti + in2, tts + 2714. ..), (CI) 



where (ao) is finite. This parametrization slightly differs from (|3.8|) , in that ai and tTj are 
defined without a factor 1/2; notice also that we do not neglect fluctuations in the direction 
of the condensate. Upon substituting (|ClD into (^TID, the functional integral becomes 



where 



Z = Jvao Vai Vtti 6{al + al + - N) e'^ (C2) 
S = Jd'r dr {(d.aof + {d,a,f + (S^vr,)^) + 



fi=x,t l,m 

l/J. 2 
— {cTod^ai - aidi_,aon2md^7!-2m.+l - T^2m+ldmuT<'2m) , (C3) 

-^^ 9iJ. 

where the indices / and m run from 1 to 2A^ — 2 and from 1 to — 1. The values of the 
couplings are the same as in ( |2.3| ). As in the body of the paper, we focus on the situation 
near the zero-temperature phase transition, which generally occurs at some g = gdl^i)- We 
also assume that the anisotropy is small, and make all computations to leading order in 7 
only. 

The large A^ expansion proceeds along the same lines as for square-lattice antiferromag- 



nets We introduce the condensate value, a, via 

ao = VN(r + ao (C4) 

and impose the constraint by introducing the Lagrange multiplier A into the functional 
integral. At A^ = cxd, the saddle-point equation gives 

2 9c ~ 9 /rfr-\ 

a = , (C5) 

9c 

where Qc is the same as in the isotropic case. 

We first consider how the ratio of the stiffnesses scales with g — g^. For this we observe 
that the only anisotropic piece in ( |C3|) which survives at infinite A^ is {■j^/2g^) cr^ (5^cri)^. 
Hence, at N = 00, the propagators for n fields remain the same as in the isotropic 0{2N) 
model 

Gtt = ^ — 5~ 27^' 

2 uji + 4k^ 
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while the propagator for the ai field acquires a correction linear in 7^ 



where cq = \J p\/x\- Eqns (|C^ ) and (|C7| ) identify (upto an overall factor) and G^-^ as 
the transverse and longitudinal propagators of gapless spin-wave excitations in the ordered 
state. Each of the propagators can now be reexpressed in terms of the fully renormalized 
spin-stiffness and spin susceptibility by collecting 7- independent terms, which are the same 
for both propagators. Comparing then the two expressions, we obtain 

= 2^^ = ,v (C8) 

where p and x ci-re now fully renormalized T = spin-stiffness and spin susceptibility respec- 
tively. We see from that while the relative difference of the bare stiffnesses is 7^, the 
ratio of the renormalized stiffnesses contains the extra factor cr^ and therefore tends to zero 
as the system approaches quantum phase transition point. We now introduce the crossover 
exponents 0i and 02 in the same way as in Sec |I|. Namely, we decompose 7^ into their trace 
and traceless parts as 

7x = 7i + 72 7r = 7i - 272 (C9) 

and define 0i and 02 as 



^^ = 7i(o)-^^ + 72(e 

X\\ - XL 



J) 



71 (O)-"^^ -272(0)-"^ (CIO) 



where O is the Josephson correlation length measured in lattice units. As g approaches 
Qc, O behaves as ~ (1 — Qx/gc) '^- At = cxo, z/ = 1. Using (|C5|), we then obtain 

01 = 02 = 1- 

Our next step will be to calculate the 1/A^ corrections to the crossover exponents. The 
corresponding diagrams are presented in Fig|^. It is not difficult to show that the polarization 
operator at T = 0, Il*{k,iuj), has no 7-dependent corrections to the leading order in 1/A^ 
and we therefore can use the same expression as in the isotropic case 



n*(fc, H = , f + - {af (Cll) 



We will also need the expression for the propagator for the fiuctuating component of the 
field along the direction of the condensate: 

^ T ujl + clk^' {^~~g ^""^^ Ii*{k,iu) ul + clk^ ' ^^^^^ 
and expressions for (a)^ and S,j with logarithmic accuracy to order 1/A^: 
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9c- g 

. 9c , 
'9c -9' 
. 9c , 



1 + ^TTT log ■ 



16 



log 



9c- 9, 
9c 



9c- 9, 



(C13) 



The evaluation of the diagrams is now straightforward. Collecting the contributions from 
all diagrams on Fig. |^ and restricting to only logarithmic contributions to order we 
obtain after some algebra 



Pll - P± ^ 9c- 9 

P± 9c 
XII -X± _ 9c-9 



Ix [1 + 
lr(l + 



GAL 

16L 
5^ 



X± 9c 

where L = log(l — 9/ 9c)- We then use ( |C9| ) and find 

112 



+ 7r 
+ Ix - 



16L 

ISTT^iV 

32L 



157r2A^ 



(CM) 



1 + 



1 + 



32 



37r2A^ 



(C15) 



These values for the crossover exponents coincide with Eqn (p.8|) obtained by other means 
in Sec. |I| 

Our next move will be to compute, in the 1/A^ expansion, the critical exponent for the 
order parameter. Using the definition ( p..l2| ) and Eqns ( |C1| ) and ( |C4| ) , we obtain 



Nn = SZs 



N 



(C16) 



where 0{1/N) stands for regular 1/N corrections from the other components of 2;-field which 
can be neglected in the calculations of the critical exponents. Using then ( pi2| ) we find 



iuj) ujI + clk'^ I 



(C17) 



Performing the integration with the logarithmic accuracy and exponentiating the result, we 
obtain 



where 



No = SZsiaY 



e = l + 



(C18) 



(C19) 



Expressing now (a) in a conventional way as (cr)^ = (Ag/gc)^, where /3 = 1 — A/Ntt^, we 
find Nq ~ (Ag/gc)^, where 



(C20) 
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We will also need the result for the staggered spin susceptibility at T = 0. Using 
1.11| ), (|C1|) , ( p^) and the result obtained in Section |3 that the number of transverse spin- 



wave modes in the ordered state is A^^^^, = 2iV, we find that the transverse spin susceptibility 
at = oo and k <S is related to the propagator of the z-field 

where index b labels transverse spin components. The relevant corrections are now the 
same as in coUinear antiferromagnets. Using the results of Ref p5|, we find that they can 
be completely absorbed into the renormalization of A^^o and p±. Assuming, as in the bulk of 
the paper, that the anisotropic term in the action will transform p± and c± into ps and c 
given by (|4.6|), we obtain 

TV^ 1 

Finally, we calculate in 1/A^ expansion, the critical exponent f]. For this, we consider the 
behavior of dynamical spin susceptibility at T = right at the transition point, g = Qc- At 
this point (a) = 0, and the spin-spin correlation function is again related to the polarization 
operator of z— fields. At A^ = cxo and g = Qc, the polarization operator is given by the 
first term in (|C11| ). To first order in 1/A^, we have to consider both self-energy and vertex 
corrections in the polarization bubble (Fig §). They both contain logarithms of external 
momentum. Evaluating the diagrams in Fig ^ with logarithmical accuracy, we obtain: 



where An = I + 0{l/N), and 
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FIGURES 



FIG. 1. a. Self-energy corrections to spinon propagator. The heavy sohd hne is a full spinon 
propagator, and the wavy line is an inverse polarization operator for spinons. The analytical 
expression for the vertex is given by (|3T^ ). b. Dia grams which contribute to current-current cor- 



relation functions < K^K^ > and < J^J^ > to first order in and to first order in 7^. The side 
vertices (shaded) in the bubbles are 2T"-^ kfj_{l + j^) / for U{1) response and 2(5"^/c^(l-|-7^/A^)/5f^ 
for SU{N) response. 



FIG. 2. Diagrams which contribute to renormalization of the F vertex (3.17) to order 1/A^. 



FIG. 3. Diagrams which contribute to the renormalization of the polarization operator to first 
order in 1/A^ at 7^ = 0. As in Fig.|^, heavy solid lines are full spinon propagators. 

FIG. 4. Second-order self-energy corrections to magnon propagators due to cubic vertices. 
Notice that cubic terms always produce anomalous self-energy terms. 

FIG. 5. Self-energy diagrams to order 1/A^ for the propagator of ai field in the ordered state 
at T = 0. Solid line is the propagator of ai given by ( |C7| ), dashed line represents the condensate 
a, wavy line is inverse polarization operator, and heavy solid line is the propagator of (Tq field 
introduced in (^: G^, = G^{1 - {2/gr)(J^G^/U*), where G^ is given by (^). 
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